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Abstract 

The control of aerial gymnastic maneuvers is challenging because these maneuvers 

frequently involve complex rotational motion and because the performer has limited 

control of the maneuver during flight. A performer can influence a manuever using a 

sequence of limb movements during flight. However, the same sequence may not 

produce reliable performances in the presence of off-nominal conditions. How do 

people compensate for variations in performance to reliably produce aerial 

maneuvers? In this report I explore the role that passive dynamic stability may play 

in making the performance of aerial maneuvers simple and reliable. 

I present a control strategy comprised of active and passive components for 

performing robot front somersaults in the laboratory. I show that passive dynamics 

can neutrally stabilize the layout somersault which involves an "inherently unstable" 

rotation about the intermediate principal axis. And I show that a strategy that uses 

open loop joint torques plus passive dynamics leads to more reliable I 1/2 twisting 

front somersaults in simulation than a strategy that uses prescribed limb motion. 

Results are presented from laboratory experiments on gymnastic robots, from 
dynamic simulation of humans and robots, and from linear stability analyses of 

these systems. 
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Abstract 

The control of aerial gymnastic maneuvers is challenging because these maneuvers 
frequently involve complex rotational motion and because the performer has limited 
control of the maneuver during flight. While a performer can execute a sequence of 
limb movements during flight to produce a desired maneuver, the same sequence may 
not produce reliable performances in the presence of off-nominal conditions. How 
do people compensate for variations in performance conditions to reliably produce 
aerial maneuvers? It is possible that people sense errors in performance and actively 
compute appropriate responses to them to produce reliable maneuvers. However, 
it is also possible that the maneuver is inherently stable, that the body naturally 
compensates for variations, and that the athlete does little active computation. 

In this thesis I explore the role that passive dynamic stability may play in making 
the performance of aerial maneuvers simple and reliable. I consider the control of 
the tucked somersault, the layout somersault, and the 11/2 twisting front somer- 
sault. I present a control strategy comprised of active and passive components for 
performing robot front somersaults in the laboratory. I show that passive dynamics 
can neutrally stabilize the layout somersault which involves an "inherently unstable" 
rotation about the intermediate principal axis. And I show that a strategy that uses 
open loop joint torques plus passive dynamics leads to more reliable 11/2 twisting 
front somersaults in simulation than a strategy that uses prescribed limb motion. Re- 
sults are presented from laboratory experiments on gymnastic robots, from dynamic 
simulation of humans and robots, and from linear stability analyses of these systems. 
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Chapter 1 



Introduction 



In both the 1992 and 1994 winter Olympics, the final jump performed by the men's 
champion freestyle skier was a quadruple twisting, triple somersault. In the per- 
formance of these 'jumps' the athletes skied off of a ramp from which they soared 
approximately 45 feet into the air, remained aloft for approximately 3.0 sec, rotated 
three times about a horizontal axis and four times about their body vertical axis, and 
finished by landing squarely on their feet so that they could continue skiing down the 
hill in a controlled fashion. It is incredible that these athletes can perform a maneuver 
like this with such accuracy. How do they do it? 

Two related issues that make these performances challenging are the controllabil- 
ity and stability of aerial maneuvers. Controllability refers to the ability of an athlete 
to influence the outcome of an aerial maneuver once it has been initiated. An airborne 
performer can not apply any external forces or torques to the body. So how can a 
performer control his body orientation? Previous research of this topic has revealed 
movement techniques that performers can use to produce complicated aerial maneu- 
vers. These techniques involve moving the limbs to reconfigure the body during flight 
(Figure 1-1). However, while a prescribed set of body configurations can produce a 
desired aerial maneuver, this sequence may not lead to reliable performances. 

The stability of aerial maneuvers concerns their reliable performance in the pres- 
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• Change inertia 
about an axis 





• Zero momentum 
reorientation 




• Reorient principal axes 




From Frohlich 1980. 



Figure 1-1: A limited ability to influence the outcome of a ballistic maneuver arises 
from the relative movement of the limbs and torso during flight. One way to influence 
the outcome is to change the moment of inertia about an axis in order to change the 
rotation rate about that axis. The standard being the ice-skater's spin. A second way 
to influence the outcome is to use momentum- free rotations [Smith 67, Frohlich 79]. 
This technique allows a structure to be reoriented while maintaining zero angular mo- 
mentum by performing a sequence of limb movements. A third way is to reconfigure 
the system so that the principle axes of inertia are reoriented relative to the inertially 
fixed angular momentum vector. This allows sharing of momentum between principal 
axes. This procedure is frequently used to introduce or remove twist in somersaults 
[Batterman 68, Frohlich 79, Yeadon 84]. The existence of these mechanisms makes 
it possible to actively adjust the outcome of an aerial maneuver once it has been 
initiated. 
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ence of off-nominal conditions. Stability is important because the accumulation of 
even small errors during the relatively long flight times of an aerial maneuver could 
result in a disastrous outcome. Making a performance reliable requires that the ath- 
lete compensate for inaccuracies in movement, variations in equipment, and external 
disturbances. It is possible that people produce reliable maneuvers by sensing these 
variations and actively computing appropriate responses to compensate for them. 
However, the complexity of this feedback control approach would appear to place 
great demands on the athlete's perceptive and motor control abilities. Is there a 
more simple approach to producing reliable maneuvers that does not require active 
compensation by the athlete? 

The focus of this thesis is on the use of passive dynamic stability as an alternative 
or a complement to active control for producing reliable aerial maneuvers. Passive 
dynamic stability means that a maneuver is inherently stable by virtue of the natural 
dynamic interaction of the limbs and body. The precise limb movement that a per- 
former makes during a maneuver will depend upon not only his motor activity but 
also on environmental forces. Is it possible that the passive forces that arise due to off- 
nominal conditions could provide a built-in correction? If so, passive behavior could 
automatically compensate for errors in initial conditions, in control movements, or 
from external disturbances. If this were the case then the demands for active control 
by the athlete could be dramatically reduced. A passive control strategy is appealing 
because it could relieve the human performer of sensing small variations in movement 
and computing control responses fast enough to produce accurate, stable maneuvers. 
Instead, the athlete may need only to "play back" a pre-recorded set of motor actions. 
This feed forward command combined with the passive dynamic response of his body 
may allow the maneuver to "unfold" on its own. 

To see if passive dynamic stability could play a role in aerial performances, I 
consider the control of three gymnastic maneuvers, the tucked front somersault, the 
back layout somersault, and the front somersault with one and a half twists. I show 
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that passive dynamics can play a significant role in making these maneuvers reliable. 
Neutrally stable passive dynamics make the control of the tucked somersault simple 
since an active controller need only consider rotation about the major principal axis. 
The layout or straight body somersault is considered to be inherently unstable because 
it requires rotation about the middle principal axis of inertia, an unstable equilibrium 
for a rigid body. I show that the layout somersault can passively be made neutrally 
stable by tuning the compliance of the shoulders of the performer. Finally, I present 
results that suggest that a passive, compliant model of the human body that uses 
open loop torque control can produce more reliable 11/2 twisting somersaults than 
a model that uses prescribed limb motion. I obtained these results using simple 
analytic models, non-linear dynamic simulation, and laboratory robots. Next, I briefly 
describe the results of experiments on each of these maneuvers. 

The somersault is a maneuver in which a performer jumps into the air and rotates 
once about a side-to-side axis before landing on the ground. The main requirement for 
the somersault is to land in a balanced manner. This in turn requires the performer to 
land with a precise body attitude. The tucked somersault exhibits passive directional 
stability in rotation since it involves rotation about the major principal axis of inertia. 
This means that imprecise initiation of a somersault will not dramatically affect the 
orientation of the spin axis during flight. However, avoiding over-rotation or under- 
rotation of a somersault about the spin axis may require compensation of rotation 
rate during flight. Active control of the inertia can be used to correct errors in 
somersault rotation rate. I present results from somersault experiments using a 60 
lb, one meter tall laboratory robot that runs on two springy legs (Figure 1-2). The 
robot was programmed to initiate the somersaults using a pre-programmed pattern 
of action. To avoid large tilt and twist angles of the 3D Biped during take-off and 
landing we use a wide double stance of the robot and insure that the feet touch down 
simultaneously. During flight the robot actively controls rotation rate by "tucking" 
or "untucking" its legs to manipulate the robots inertia. The robot actively adjusts 
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Figure 1-2: Photograph of the 3D Biped Robot. 

the position of its feet relative to its center of mass just prior to landing. This element 
was necessary to compensate for errors in the estimated landing time of the robot. 
On its best day the robot did successful somersaults and continued running on 7 out 
of 10 attempts. 

The layout or straight body somersault is considered to be inherently unstable 
because it requires rotation about the middle principal axis of inertia, an unstable 
equilibrium for a rigid body. A rigid body that is somersaulting about the middle 
principal axis will always exhibit a sequence of half twists about the minor principal 
axis. Despite this fact, athletes regularly perform this maneuver with apparent ease. 
Previously, biomechanics researchers have assumed that the athlete senses the insta- 
bility of the maneuver and actively compensates for it with movements of the arms 
and body [Nigg 74, Hinrichs 78, Yeadon 90]. 

I show that the layout somersault can be a passive, neutrally stable maneuver. 



CHAPTER 1. INTRODUCTION 17 




Figure 1-3: Photograph of the mechanical doll used for experiments with layout 
somersaults. The doll could consistently perform stable, triple layout somersaults. 

Passive stabilization of the layout somersault results from the natural dynamic inter- 
action of the limbs and body during movement. Stabilization arises from the inherent 
tendency of the arms to tilt in response to twisting movement of the body. The arm 
tilt forces the principal axes of the system to move in a direction that compensates 
for tilt and twist errors. This built-in correction eliminates the divergent tendency of 
the system as long as the compliance of the shoulders cancels the unstable centrifu- 
gal forces on the arms. I verify this result with linear stability analysis, non-linear 
dynamic simulation, and experiments on a human-like doll built and tested in the 
laboratory (Figure 1-3.) The doll has spring-driven arms but has no other control 
system, sensors, or actuators. The doll routinely exhibits triple somersaults about its 
middle principal axis without twisting. 

While twisting is to be avoided in the layout somersault, it is a feature in other 
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Figure 1-4: Images arranged in right-to-left, top-to-bottom order from a dynamic 
simulation of a 1 1/2 twisting front somersault. We found that open loop torque 
control of a compliant model led to relatively reliable performances. 



tricks. I present results from dynamic simulation of a human performing a front 
somersault with one and a half twists. The twisting maneuver is started from a front 
somersault by reconfiguring the body mid-maneuver. While a sequence of prescribed 
limb movements can be found to produce a twisting maneuver, simulation results 
show that the performance is sensitive to small variations in initial conditions. If, 
instead, open loop torque control is used with a compliant, passive dynamic model of 
the human the maneuver reliability can be significantly improved. 

I also present results from dynamic simulation of a front somersault with one half 
twist performed by the 3D Biped. We found that to make the simulated 3D Biped 
perform a satisfactory front somersault with one half twist we had to add weight to 
the robot to make its inertia more like that of a human. 

The ultimate goal of this research is to develop a theory of passive dynamic, 
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aerial maneuvers. By showing that passive dynamics can play a significant role in 
the performance of reliable aerial maneuvers we argue that they reduce the need for 
active control from the athlete. The performance of a maneuver may be simplified 
by relegating some responsibility for control to the natural, mechanical behavior of 
the body. The appropriate use of passive dynamics may have the added benefit 
of producing natural looking, coordinated movement. Perhaps athletes and other 
people use a performance strategy that seeks to maximize passive dynamic behavior. 
It is difficult to know what control strategies people may or may not use, but in 
the laboratory we can examine the feasibility of a strategy by testing it in a real or 
simulated system. 

1.1 Background 

Relevant background material for the study of gymnastic maneuvers comes from 
fields such as biomechanics, biology, robotics, aeronautics, and astronautics. Several 
researchers have explicitly studied the performance of gymnastic maneuvers. These 
studies have revealed the salient features of known gymnastic techniques for produc- 
ing maneuvers. Some robotics researchers have studied gymnastic maneuvers using 
dynamic simulation and/or laboratory robots in order to develop strategies for con- 
trol. The study of passive dynamic stability has roots in the fields of aeronautics 
and astronautics. It is common for airplanes and spacecraft to be designed to exhibit 
passive dynamic stability. The study of human locomotion provides an inspirational 
example of passive dynamic stability and the rich behavior that it can produce. In 
the following sections, I briefly discuss background material from each of these fields. 

1.1.1 Gymnastic Maneuvers 

Most researchers investigating the control of aerial maneuvers have been concerned 
with explaining the physics of the maneuvers. During the aerial phase of a maneuver 
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the performer can not produce any external forces on the system so the trajectory of 
the center of mass and the angular momentum are fixed from the moment of take- 
off to the moment of landing. Knowing this, early researchers were compelled to 
find explanations of two aerial maneuvers that at hrst appeared to violate conser- 
vation of angular momentum: 1) a cat when dropped with no net angular momen- 
tum will right itself before landing on the ground [Marey, Kane 69] and 2) spring 
board divers who leave the board with rotation only about a side-to-side (somer- 
sault) body axis can subsequently initiate rotation about their head-to-toe (twist) 
body axis [Batterman 68, Frohlich 79, Yeadon 84]. Where did the extra angular mo- 
mentum come from? Resolving these apparent discrepancies led researchers to find 
movement techniques that could be used to perform these and other interesting ma- 
neuvers. These techniques involved reconfiguring the body in flight (Figure 1-1). 

Takashima [Takashima 90] studied high bar maneuvers. He used the control of 
rotation rate about the somersault axis to produce a balanced landing of a simulated 
human dismounting from the horizontal bar. His algorithm could produce accurate 
tucked, multiple-somersault dismounts. In flight, he used the feedback control of 
posture to control rotation rate. He used a combination of feed forward and feed 
back control to execute the landing. 

While they appear to be similar maneuvers, the tucked somersault and the layout 
somersault are dynamically very different. The tucked somersault involves rotation 
about the major principal axis of inertia, a stable mode of rigid body rotation. Rota- 
tion about the major and minor principal axes of inertia exhibits a limited form of sta- 
bility called directional stability. Directional stability refers to the fact that the spin 
axis will maintain roughly the same inertial orientation when deviated slightly from 
that orientation. The layout somersault involves rotation about the intermediate prin- 
cipal axis of inertia, an unstable mode of rigid body rotation [Crandall 68, Hughes 86]. 
When a rigid body is spun about its intermediate axis, that body axis will exhibit 
large excursions from its initial orientation. These cyclic excursions are a series of 
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nearly 180° rotations away Irom and returning to the initial axis orientation. The 
different modes of rotation can be easily demonstrated to the reader by spinning a 
video tape or a book (with a rubber band around it) about each of its axes of sym- 
metry. You will see that when the body is spun about the major and minor principal 
axes the spin axis roughly maintains its inertial orientation. However, when the body 
is spun about its middle principal axes it will exhibit a sequence of twists about the 
long body axis. This rigid body instability has lead researchers to conclude that the 
layout somersault is inherently unstable. 

Nigg recognized that the layout somersault may be unstable since it involved 
rotation about the middle principal axis [Nigg 74]. Using cinematographic techniques 
Hinrichs [Hinrichs 78] measured the body configuration of an athlete performing a 
layout somersault. He confirmed that the layout somersault involved rotation about 
the middle principal axis of inertia. He hypothesized that the athlete made small 
corrective movements of the arms and torso in flight to stabilize the somersault. 

Yeadon [Yeadon 84] used a combination of cinematography and dynamic simula- 
tion to study the control of aerial maneuvers. In his research he developed a mass 
properties model of the human form that provided an estimate of inertial parameters 
from anthropometric measurements. He also filmed highly skilled athletes performing 
complex aerial maneuvers and digitized this data to determine the body attitude, con- 
figuration, and angular momentum during flight. Then he used the inertia parameters 
and digitized configuration and momentum data as input to a dynamic simulation of 
the human body during flight. This dynamic model computed the gross body atti- 
tude during flight as a function of the measured internal configuration and angular 
momentum. Using this system, Yeadon could numerically study the effect of changes 
in body configuration on the performance of complex aerial maneuvers. 

Yeadon found that by piking or arching (bending at the waist in the sagittal 
plane) during a layout somersault an athlete can change his or her inertia enough 
to make the somersault axis an axis of maximum inertia. This in turn implies that 
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arched /piked somersaults are passively stable. Because many athletes perform layout 
somersaults with arch, it is possible that this is done to stabilize the maneuver. 

However, in Yeadon's analysis of a double layout somersault performed by Carl 
Furrer, the 1982 World Trampoline Champion, he found that Furrer was in fact ro- 
tating about the middle principal axis of inertia during nearly all of the maneuver. 
Furthermore, dynamic simulation of this layout somersault exhibited the character- 
istic twist instability of rotation about the middle principal axis while the actual 
human performance exhibited no such instability. This difference implied that the 
small digitization errors in translating him configuration data to the simulation were 
responsible for the change in performance, a fact that would point to an inherently 
unstable system. These results suggest that the athlete uses some form of stabilization 
to perform the layout somersault. 

Yeadon proposed a specific technique for stabilizing the layout somersault. He 
realized that asymmetric movement of the arms in the frontal plane could be used to 
stabilize rotation about the middle principal axis of inertia. He designed a stabilizing 
feedback controller that used the sensed twist angle as the feedback signal to drive 
the arm abduction/adduction angular rates. Using a linear model, Yeadon found 
the athlete would have to respond to a growing twist instability within 0.28 of a 
somersault (ss 200ms) [Yeadon] in order to maintain stability. 

The ballistic nature of many gymnastic performances makes an open-loop (feed 
forward) component a likely part of any control strategy. An open loop strategy is 
one in which the performer's control motions are not derived from the current state 
of motion but are instead produced from a pre-programmed pattern of action. Since 
important parameters of ballistic motion are fixed from takeoff, an open-loop strategy 
is required to anticipate the maneuver in order to set up appropriate initial conditions. 
Raibert and Hodgins programmed a planar biped robot to perform front somersaults 
in the laboratory using a combination of open loop and feedback control strategies. 
The biped robot used an open loop strategy for initiation of the somersault and the 
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majority of flight. A feedback controller was used to accurately place the feet just 
prior to touch down. Using this combination of feed forward and feedback control, the 
robot could successfully perform the front somersault followed by balanced running 
on 90% of the trials. 

1.1.2 Aeronautics, Astronautics and Celestial Mechanics 

Some modern aircraft are designed to require active control for stabilization. These 
aircraft use digital computers to manipulate the aircraft control surfaces to render the 
planes flyable. This inherent instability is tolerated to provide highly maneuverable 
aircraft. However, passive dynamic stability is commonly built into general aviation 
aircraft and spacecraft. More precisely, they have an equilibrium condition like flying 
straight and level for an aircraft, or spinning about the major principal axis for a 
satellite that is stable so that the craft can tolerate external disturbances without 
diverging from the stable equilibrium condition. The origin of this topic has its roots 
in celestial mechanics. 

The moon always presents the same face to the Earth. However, it does not do 
so exactly. Librational stability refers to the stability of the oscillations of the moon 
about its center of mass as it circles the Earth. Galileo was the first to notice these 
oscillations. Newton conjectured that a reason for this behavior would be that the 
moon was elongated towards the Earth. However, it took Louis Lagrange to develop a 
mathematical theory describing this phenomenon [Lagrange]. Using a linear analysis, 
Lagrange derived a set of four inequalities involving the inertia of the moon that must 
be satisfied for the moon to exhibit stable librational motion. 

In 1885, Henri Poincare [Poincare] realized that a linearized analysis could not 
be conclusive in determining librational stability. This result gave rise to the use of 
the Hamiltonian as a Lyapunov function candidate in determining Lyapunov stability. 
This approach established the stability of a satellite configuration called the Lagrange 
satellite but it could not establish the stability of the Delp satellite which was also 
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stable according to linear analysis. 

The only other example that the author has found of stable rotation about the 
intermediate axis is from a study of large cellestial structures. Duncan and Levison 
[Duncan 89] simulated the behavior of a self-gravitating system of 2048 bodies in 
order to determine if it was stable. They found an example of a simple spherical 
system initially in dynamic equilibrium that experienced an instability producing a 
final equilibrium state of stable rotation around the intermediate principal axis of 
inertia of the system of particles. This result was considered noteworthy because it 
conflicted with the rigid body analogy of unstable rotation about the middle principal 
axis. No explanation of the source of stability was offered. 

1.1.3 Ballistic Walking 

While the significance of passive dynamic stability is recognized in studies of celestial 
structures and in the design of aircraft and spacecraft its relevance to the control of 
movement in biology and robotics is only beginning to be explored. One possible 
reason is that these former examples typically involve the stability of an equilibrium 
configuration, i.e. no accelerations. Animals, people, and robots frequently move 
with significant accelerations. The stability analysis of non-equilibrium movement is 
a much more difficult process. Some progress has been made with the analysis of 
walking. 

Mochon and McMahon and later McGeer showed that passive dynamic stability 
may be important to human locomotion. First proposed by Mochon and McMa- 
hon [Mochon 80], a ballistic walker uses only gravity and the dynamic interaction 
of the swing and stance legs to produce a repetitive walking pattern. The passive 
pattern accounted for the folding and unfolding of the legs and the positioning of 
the foot forward. McGeer [McGeer 89] showed the viability of ballistic walking by 
building passive, planar, anthropomorphic linkages with no sensor or actuators that 
demonstrated stable walking down an incline. 
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1.2 Organization of Thesis 

The remainder of this thesis is organized as follows. 

Ch. 2 This chapter provides a review of rigid body rotation. Rigid body rotation 
provides an important simplified model for the analysis of multi-body rotational 
systems. Analysis and tools introduced here will be used throughout the thesis. 

Ch. 3 In this chapter I discuss somersault experiments with a 3D biped robot. The 
robot somersault axis is coincident with the major principal axis of inertia so the 
maneuver exhibits some passive stability properties. The robot actively controls 
landing attitude by retracting or extending its legs during the maneuver to 
change somersault rate. The robot has successfully performed front somersaults 
in the laboratory. 

Ch. 4 In this chapter I present an analysis of the layout somersault. I show that the 
layout somersault, involving rotation about the middle principal axis of inertia, 
can be passively stabilized by tuning parameters of a passive dynamic model of 
the human body. Using the simplest possible model of the layout somersault, I 
explain the fundamental dynamics of passive stabilization. 

Ch. 5 I describe layout somersault experiments with a human-like doll. These exper- 
iments verify that the layout somersault can be consistently stabilized for at 
least three and one half somersaults. 

Ch. 6 I discuss non-linear dynamic simulation of a I 1/2 twisting front somersault. 
Simulation results suggest that maneuvers using prescribed limb motion will 
not be reliable to off-nominal initial conditions. If instead, a compliant passive 
model is used with feed forward torques the maneuver can be made more reli- 
able. I also describe twisting somersault experiments with the 3D Biped robot. 
Using simulation we found that it was important that the robot have an inertia 
tensor more like that of a human to produce a human-like front somersault with 
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1/2 twist. Laboratory experiments with 3D Biped robot twisting somersaults 
were not successlul at least in part due to insufficient actuator power. 

Ch. 7 Here I summarize the results of the thesis and discuss future work. 

Appx. In the Appendix I provide Mathematica code for deriving non-linear equations 
of motion of a human model, and for analytically linearizing this model. I 
provide a definition of the parameters used in the linear model. I also provide 
a derivation of the analytic solution for rigid body rotation. 



Chapter 2 



Rigid Body Rotation 



2.1 Introduction 

Humans are multi-body systems and gymnastic maneuvers involve multi-body rota- 
tions. We can better understand multi-body rotation by understanding the dynamics 
and stability of rigid body rotation. In this chapter, I briefly discuss some properties 
of rigid body rotation. Rigid body rotational stability about a principal axis de- 
pends only on the relative magnitude of the principal inertias. This result provides a 
valuable reference for multi-body rotational stability. A linear analysis of rigid body 
rotation provides simple stability results, and identifies important non-dimensional 
parameters that will be useful in multi-body analysis. Visualization tools for rigid 
body rotation will also prove useful in understanding how to control rotational motion 
in multi-body systems. Material for this chapter is based on the text of [Hughes 86]. 

2.2 Nonlinear Dynamic Model 

The rotational motion of a free rigid body about its center of mass is governed by 
Euler's equations. 



27 
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ItUt - (I 2 - I 3 )U 2 L0 3 = 

I2LO2 - (h - h)^3^i = (2.1) 

I 3 u 3 - (7i - I 2 )loilo 2 = 

where Ij refers to the j th principal moment of inertia and uj 3 refers to the angular 
velocity about that principal axis. 

The orientation of the rigid body with respect to an inertial coordinate frame is 
described by a set of three Euler angles. The Euler angles are used to define the 
3-by-3 matrix of direction cosines, Cu, that transforms a vector described in inertial 
coordinates into one described in the body fixed coordinate system. (For the rigid 
body analyses of this chapter I assume that the principal axis frame is coincident 
with this body axis frame.) I use a '2-1-3' Euler angle sequence for this purpose. 
I borrow the names for the three Euler angles, somersault, tilt, and twist, from 
Yeadon [Yeadon 84]. I use the letters, $,0,^ to refer to the somersault, tilt, and 
twist angles respectively. 

To describe the attitude of the body in inertial space, a coordinate system initially 
parallel to the inertial reference frame is first rotated through the somersault angle 
about the inertially fixed '2' axis, then rotated through the tilt angle about the 
intermediate '1' body axis, and finally rotated through the twist angle about the 
body fixed '3' axis. Figure 2-1 depicts the $, 0, *& Euler angle sequence applied to a 
rigid human form. The rotation matrix is given by 



Ch- 



ClJrO<]> + »J>i])aJ>0aJ><I> b\$C<~) — (y\$b§ ~\~ b If b Q G <J> 



(2.2) 



The subscript hi refers to the fact that this rotation matrix rotates a vector from the 
inertial system into the frody axis system. S and C are the sine and cosine of the 
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Figure 2-1: Illustration of the somersault ($), tilt (0), and twist (\P) Euler angle 
sequence used to define body attitude. 

respective angle. This Euler angle description of body attitude has a singularity, as do 
all three parameter descriptions of body attitude. The singularity for this particular 
sequence occurs at a tilt angle of ±7r/2. The kinematic equations governing the 
evolution of the Euler angles are given by 







SxS/SqGq Gx^oqGq 1 

G\$> —S\s> 

Ox^Gq Gx^iGq 



u 3 



(2.3) 



The inverse of the above relationship is also useful. It is given by 



u 3 



G\$> Sx^Gq 

— iS\]j Gx^iGq 
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(2.4) 
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The dynamic equations of motion (2.1) are coupled so that general rotational 
motion involves all three degrees of freedom. However, three simple solutions exist to 
these equations. If a rigid body is spinning perfectly around a principal axis such that 
two of the three angular velocities are zero then the non-linear coupling terms vanish. 
The body will continue to spin about that axis without coupling to the other degrees 
of freedom. An analysis of these spin solutions will reveal that for a tri-inertial body, 
a body with three different principal inertias, only two of the spin solutions are stable 
while the third is unstable. Here stability means if the spin axis of the system is 
moved away from the nominal solution the attitude of the spin axis will not diverge 
radically from its initial orientation. The linear stability analysis provided in the next 
section reveals the required conditions for stable rotation about a principal axis. 

2.3 Linear Analysis 

In this section we linearize Euler's equations for the rigid body in order to study the 
stability of steady rotation about a principal axis. The linear equations govern the 
motion of the body relative to a reference frame that is rotating steadily about the 
principal axis of inertia. The resulting equations will be valid for small deviations of 
the system from the reference frame. This linearization process will be used again 
when we linearize multi-body rotation about an equilibrium spin condition. The 
linearizing condition will be steady somersaulting rotation about the '2' axis with 
angular velocity $. The angular velocity of the steadily rotating reference coordinate 
system is u ra . 

u ra = $z r2 (2.5) 

where i r2 is the unit vector along the '2' axis in the reference coordinate system. The 
angular velocity of the principal axis system, a3 pa , is comprised of the angular velocity 



CHAPTER 2. RIGID BODY ROTATION 



31 



relative to the reference coordinate frame, uv, and uf ra 



UJpa UJ pT ~r Lo ra 



(2.6) 



Let the attitude ol the principal axis system relative to the reference frame be de- 
scribed by the 2-1-3 Euler angle sequence of <f),6,ip. (The lower case notation indi- 
cates that these are linearized states. These Euler angles describe the deviations of 
the body from the rotating reference frame.) We can use Equations 2.2 and 2.4 to 
express the components of u pa in terms of </>, 6, ip, their derivatives and $. First using 
Equation 2.4 

C %i) S^Ce ip 

-SV C\ p C e (2.7) 

1 -S e 

where u pr denotes a column whose components are the elements of the vector u pr 
when described in the principal axis frame. To form the sum in Equation 2.6, we 
need to express all the vector components in the same frame. To do this, using 
Equation 2.2 as a model, form the rotation matrix C pr that rotates a vector from the 
reference frame into the principal coordinate frame. Then solve for io pa} 



u, 



pr 



UJpa LUp r ~\- ^jp r LU ra 



When C pr and u pr are linearized for small <f>,ip,6 we get 



u, 



pa 
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(2.8) 



CHAPTER 2. RIGID BODY ROTATION 



32 



This expression can be differentiated to find the Euler angle expression for angular 
acceleration. 



uJpa 



9 




V> 
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+ $ 





J . 




-9 



(2.9) 



Now, substituting Equations 2.8 and 2.9 into Equation 2.1 and eliminating terms that 
are second order in 6, <f>, and ip results in the following equations 



he-§{h-h-h)i>-&{h-i 2 )e = o 



h 



hil> + ${l2-h-h)6-fr{h-I 2 )il> = 



(2.10) 



These equations govern the behavior of the non-linear system in the vicinity of the 
pure spin solution. The linear states, 6, <f>, ip, describe the deviation of the body axes 
from the reference coordinate frame. 

Examination of Equations 2.10 reveals that the system is unstable in the sense 
that perturbations in <f> will result in unbounded growth in <f>. Nevertheless, a reduced 
form of stability called directional stability [Hughes 86] is possible in the subsystem 
of 6, t/j. Directional stability means that the two attitude variables 6 and ip will not 
diverge from zero if the system is perturbed slightly from the pure spin solutions. 
This means that the spin axis will continue to point in roughly the same inertial 
direction when disturbed from its equilibrium orientation. 

The equations for 6 and ip decouple from <f>. They can be written as 



Mx + $Gx + $ 2 Kx = 



(2.11) 



where 



9 
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M 



h 
h 



G 



o -(h-h-h] 
h-h-h 



K 



h-h 

h-h 



2.3.1 Conservative Gyric Systems 

To investigate the stability of the 6 — if) subsystem, I use the matrix second order 
stability theory described in [Hughes 86]. Hughes describes a second order system of 
equations of the form of Equation 2.11 as a conservative gyric system when 

M T = M > 
G T = -G 
K T = K 



The coefficient matrices (Af, G and K) are respectively associated with inertial, gyric, 
and stiffness forces. Conservative refers to the fact that the system energy is con- 
served, while gyric reflects the fact that these terms often arise in spinning systems. 
This form of equations will be present in a multi-body analysis of rotating systems. 

Hughes proves that asymptotic stability for a conservative gyric system is impos- 
sible by showing that if s is a root of the characteristic equation of 2.11 then — s must 
also be a root. Strictly left half plane poles will always have their right half plane 
counterparts. Therefore, stability, as opposed to asymptotic stability, is the strongest 
result possible for a conservative gyric system. A stable system will have all of the 
roots of its characteristic equation on the imaginary axis. 

A sufficient condition for stability of a conservative gyric system is that it be 
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statically stable, i.e. the stiffness matrix must be positive definite, K > 0. The 
stiffness matrix for the two-by-two system above is positive definite if I 2 is the major 
principal axis of inertia. That is, rotation about the major axis is directionally stable 
by virtue of its static stability. However, K > is only a sufficient condition for 
stability. Rotation about the minor principal axis is also stable. Since in this case the 
stiffness matrix is not positive definite the system is considered gyrically stabilized. 
To test for gyric stability we can check the roots of the characteristic equation of the 
second order system. The system will be considered stable if the roots are purely 
imaginary. 

These results do not preclude the possibility that asymptotic stability could be 
obtained by adding damping to a conservative gyric system. However, while damp- 
ing tends to make statically stable systems become asymptotically stable systems, 
damping also tends to destabilize gyrically stabilized systems. Hughes proves that 
statically unstable systems are destabilized if they have a positive definite damping 
matrix. 

2.3.2 Rigid Body Inertia Ratios 

Before we solve for the roots of the characteristic equation of Equation 2.11 we should 
note that we can simplify our analysis by recognizing that only the ratios of inertia 
are important to the analysis rather than the individual values of inertia. Dividing 
each equation of Equation 2.11 by the corresponding diagonal term of M results in 
the following form of M , G and K. 



M 



1 
1 



G 







1 - k x 
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K 



fa o 
o fa 



where 



, 1 2 — I 3 , 1 2 — I\ 

fa = ; , fa = 



h ' h 

A physical body can not have arbitrary values of /i,/2 and I 3 [Hughes 86]. These 
constraints are captured by the fact that 



-1< h < 1 

-1< fa < 1 

Therefore, all possible rigid body configurations can be represented on a plot of the 
parameter space of fa and fa restricted to the unit square. Figure 2-2 shows how fa 
and fa depend upon the relative size of Ji, 7 2 , and I 3 . Also included in the figure are 
schematic drawings of rectangular prisms that would have approximately the correct 
inertia ratios for selected points around the diagram. 

The characteristic equation of this system is formed with the following determinant 

det[Ms 2 + $ Gs + & K] =0 

which results in the following 



^ 4 + i,(l) 2 + i 2 = o 



(2.12) 



where 



&i = 1 + fafa 
b 2 = fafa 
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Figure 2-2: The inertia ratio diagram shows how k\ and k 3 depend upon the relative 
value of the principal inertias. k\ = 2 ~ 3 , k 3 = 2 ~ 1 . The different rectangular 
prisms located near the axes and corners of the diagram indicate an example shape 
that would correspond to the local inertia ratios. 
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For stability we require 

1. ^ = I + &1&3 > 

2. 6 2 = &1&3 > 

3. (60 2 -46 2 = (1 - ^fe) 2 > 

The only stability condition not met automatically is k\k 3 > 0. This requirement 
specifies that k\ and k 3 must have the same sign. On a plot of the {M, k3} parameter 
space stability occurs in the hrst and third quadrants. 

Figure 2-3 shows the regions of stability for the simple spin solutions of a rigid 
body. 

Equation 2.12 is simple enough that we can solve for the roots of the characteristic 
equation in closed form. Two roots are located at s = ±j$. This is the stroboscopic 
mode of rotation. This mode of rotation occurs when the body is still spinning 
perfectly about the principal axis but this axis is offset from the original orientation. 
The other two roots are located at s = ±z/, where 

v = (hh) 1 ' 2 * 

In the hrst and third quadrants of the inertia ratio diagram, v is positive and therefore 
these roots are purely oscillatory. In the second and fourth quadrants of Figure 2-2, 
v is imaginary forcing the roots of the characteristic equation to have real positive 
and negative values. In these quadrants, the unstable mode of motion is governed by 
the following equation 

x = x e s where s = (-k^f^t (2.13) 

To get an idea of how unstable the system is we compute the change in the nominal 
somersault angle, $, required for the the unstable rotational mode to grow by a factor 
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Figure 2-3: Stability diagram for the simple spin solutions of a rigid body. Rigid 
body rotational stability depends only upon the two non-dimensional inertia ratios 
k\ and k 3 . We assume the body is rotating about the principal axis corresponding to 
I 2 for this diagram. 
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Figure 2-4: All rigid bodies with inertia ratios on a single curve have an unstable 
mode that grows at the same rate. The value of a curve, indicated in the figure, is 
the somersault angle (in radians) the body must execute before the unstable twist 
mode grows by a factor of ten. The most unstable systems are those with inertia 
ratios in the corners {ki 7 k 3 } = {1,-1} or { — f , f }. This plot is symmetric about the 
origin. 



of TV. 



$ = $* 



ln(N) 



Figure 2-4 shows curves of constant $ on the {ki 7 k 3 } axes for N = 1 0. 



(2-14) 



2.4 Rotational Stability of the Rigid Human Body 



For any particular configuration of the human body, we can solve for the orientation 
and magnitude of the principal axes of inertia. This allows us to compute the cor- 
responding inertia ratios which in turn provide a valuable reference for the stability 
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Figure 2-5: Tuck, pike, and layout positions for a human. 



of the human body in rotation about the principal axes. Figure 2-5 shows from left 
to right a human in the tuck, pike and layout positions. Figure 2-6 shows the inertia 
ratios of a human performer moving from a tuck position through a pike position to 
a layout position. To make this figure, I assumed that I 2 was the principal inertia 
along the somersault axis of the human. This analysis shows that for a rigid body, 
tuck and pike somersault configurations are stable in rotation about the somersault 
axis while steady rotation in the layout position is unstable. In the pike and tuck 
positions the inertia ratios are in the stable upper right quadrant of the inertia ratio 
diagram while the inertia ratios of the layout somersault are in the unstable lower 
right quadrant. 
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Figure 2-6: Inertia ratios of a human performer for a sequence of configurations 
connecting tuck, pike and layout positions. The tuck and pike positions are stable 
(for a rigid body) while the layout is unstable. 
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2.5 A Map of Rigid Body Rotation 

The linear analysis we have used so far is restricted to three solutions for rigid body 
rotation, pure spin about each of the principal axes. A general solution for the torque- 
free motion of a rigid body tumbling in space would be useful for understanding the 
range of behavior a rigid body can exhibit between these three special solutions. 
Solving for such a solution is one of the classic problems of dynamics. The con- 
tribution that this solution offers today is a concise description of the states and 
non-dimensional parameters that govern rigid body rotation (see Appendix A. 4). In 
addition, the analytic solutions give rise to elegant geometric interpretations of rigid 
body rotation that help to provide intuition of this otherwise complex movement. In 
this section, I present one form of geometric tool derived from the analytic solution 
that is useful in visualizing rigid body rotation. I think of this tool as a map of rigid 
body rotation because it shows graphically two of the three Euler angle trajectories 
involved in non-linear rotational motion of a rigid body. This map not only makes 
clear the stable and unstable axes of rotation but also shows two distinct regions 
of qualitatively different motion. A performer can exert control over his rotational 
motion by moving within and inbetween these two regions. 

Two integrals of motion are used to define a map of rigid body rotation. During 
flight, a rotating rigid body must conserve angular momentum, h, and kinetic energy, 
T. The magnitude of both of these quantities can be written as the equation for an 
ellipsoid in angular velocity space 

h 2 = ]\u\ + I 2 u> 2 2 + T\ujI (2.15) 

2T = I x u? x + W 2 + I 3 u 2 3 . (2.16) 

Assume, with no loss of generality, that the angular momentum vector is aligned with 
the inertially fixed '2' axis. Using the direction cosine matrix, (2.2), we can solve for 
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the components of h in the principal coordinate frame. 

h = hc 1 i pl + hc 2 i p2 + hc 3 i p3 

where {ci, c 2 ,c 3 } are the direction cosines of h and i PJ is the j th unit vector in the prin- 
cipal axis system. Alternatively, the principal axes components of angular momentum 
can also be written as follows 

h = JlCUiZpi + I 2 co 2 i p2 + h^3ip3 

Thus we have established that 

hc-i = Jicui, hc 2 = I 2 to 2l hc 3 = I 3 lo 3} 

The ellipsoid describing angular momentum (2.15) can now be written as the equation 
of a sphere in direction cosine space. 

c\ + c\ + c\ = 1 

Similarly the kinetic energy can be written as the equation of an ellipsoid. 

T 1 + T 2 + T 3 = l^ 

For a hxed kinetic energy and angular momentum, the direction cosines must si- 
multaneously lie on the surface of both the momentum sphere and energy ellipsoid. 
Therefore, the intersection of these two shapes describes a trajectory in direction 
cosine space that the body must 'follow'. 

Figure 2-7 shows a sample map for the possible rotational trajectories of a 'rigid' 
human body. Each trajectory shown on the sphere corresponds to a different kinetic 
energy of rotation. The axes of the sphere are the direction cosine axes. The sphere 
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Figure 2-7: A map of rigid body rotation for a human in the layout position. The 
axes of the sphere and the principal axes of the body remain parallel as the body 
rotates. The inertially fixed angular momentum vector (black) paints trajectories 
onto the surface of the sphere as the sphere rotates. Each trajectory on the sphere 
corresponds to a different rotational energy. The trajectories indicate the tilt and 
twist angles of the body as it rotates. This map does not include somersault angle 
and does not show the time dependence of the tilt and twist Euler angles. 



rotates with the body so that its axes remain parallel to the principal axes of the 
body. The angular momentum vector is shown protruding from the sphere. The 
sphere (and body) must move so that the (inertially fixed) angular momentum vector 
remains in the 'slot' that is appropriate for the given kinetic energy. 

The stable and unstable axes of rotation are obvious from this map. The stable 
principal axes are surrounded by trajectories that enclose the axes while the unstable 
middle principal axis shows trajectories that converge then diverge from that axis. 
This map also shows that the trajectories are divided into two qualitatively different 
regions. This division is based on the rotational kinetic energy. 
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For a fixed angular momentum, the energy of rotation is bounded above and below 
by the rotational energy in pure spin about the minimum and maximum principal 

axes respectively. 

h 2 h 2 

< T < 



21 ~ ~ 21 

'- J - L vriax '- J - L vmn 

Between these extremes exists a continuum of energy levels of rotational motion. 
The dividing point among these trajectories is the energy required to spin about the 



middle principal axis, T 



2/„ 



-. Those trajectories with higher energy involve monot- 



ically increasing (or decreasing) twist while those with less energy involve oscillatory 
twisting. In Figure 2-7, the trajectories which enclose the minimum principal axis 
(head-to-toe axis) of the performer involve monotonic twist while those that enclose 
the maximum principal axis of inertia involve twist angles that oscillate between 
and ±7r depending on which region the trajectory is located in. 

We can read the twist and tilt Euler angles of a trajectory from this map in an 
intuitive way. Using Equation 2.2 we can derive an explicit relationship between the 
direction cosines and the Euler angles as follows 
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Drawing the components of the angular momentum vector as shown in Figure 2-8, we 
realize that the spherical coordinates of the unit angular momentum vector are defined 
by the tilt and twist Euler angles, {0, ^}. The tilt angle, is the (negative) latitude 
of the energy curve and the twist angle, *& , is the (negative) longitude. Therefore we 
can simply read the tilt and twist Euler angles of the body from the polar coordinates 
of the angular momentum vector on this map. For example, the attitude of the body 
shown in Figure 2-7 is approximately = 3°, \& = 0. 
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Figure 2-8: The spherical coordinates used to define orientation of the unit angular 
momentum vector in direction cosine space are the negative of the tilt and twist Euler 
angles. 

2.6 Dynamic Simulation Environment 



The analytic and experimental results presented in this thesis are complemented with 
numeric results from non-linear dynamic simulation of the creature in question. This 
section briefly describes the simulation environment used to compute the motion and 
produce graphic output from dynamic simulations. 

The simulation environment consists of three parts that allow us to compute the 
movement of a creature, make a movie of its motion using computer graphics, and 
analyze its motion using time histories of simulated data. Using software developed 
in this laboratory, the Leg Laboratory at MIT, we have integrated the simulation en- 
vironment so that the three parts are all generated automatically from a single input 
hie that describes the shape, topology, and mass properties of the system. The part 
that computes the motion does so by numerically integrating the non-linear equations 
of motion for a given creature. The non-linear equations of motion are derived us- 
ing the commercially available software package, SD-FAST. This package derives the 
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equations of motion for a rigid link, multi-body system using the description file as 
input. The equations of motion are automatically incorporated into a dynamic simu- 
lation that handles user interface to the software. The computer graphic and analysis 
software is also generated automatically from the description of the creature allow- 
ing us to animate movement from simulation data using either a simple "working" 
picture or a higher quality, computer graphic image. The only part of the simulation 
environment that is not automatically generated is the control software. This soft- 
ware is used to specify the desired behavior of the system. It does so by computing 
joint torques or other actuator inputs that we assume the creature to have. Passive 
forces like those due to springs, dampers and gravity are input to the equations of 
motion in a manner similar to active joint torques. The ability to quickly produce 
dynamic simulations of a variety of creatures has made this simulation software a 
useful research tool in the investigation of machine, human or animal movement. 

2.7 Summary 

Material in this chapter is derived from the text of [Hughes 86]. In this chapter, I 
reviewed the dynamic and kinematic equations of rigid body rotation. Linear equa- 
tions of motion describe rigid body rotation in the neighborhood of the equilibrium 
solution of pure spin about a principal axis. A stability analysis of these equations 
shows that rigid body rotation is stable about the major and minor principal axes but 
is unstable about the intermediate axis. The linear stability analysis can be simpli- 
fied by the use of the non-dimensional inertia ratios, k\ and k 3 . A stability diagram 
simply demarcates stable and unstable regions in the &i, k 3 parameter space. I show 
the values of k\ and k 3 that a (rigid) human body would assume in the tuck, pike, 
and layout positions in order to provide a stability reference for the rotating human 
body. A geometric map derived from the closed-form solution of rigid body rotation 
captures the range of possible rotational modes. This map intuitively shows the stable 
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and unstable axes of rotation. It also demarcates between two regions of qualitatively 
different rotational motion. 



Chapter 3 



Robot Tucked Somersaults 



3.1 Introduction 

The somersault is a maneuver in which a performer jumps into the air and rotates 
once about a side-to-side axis before landing on the ground. The main requirement of 
a successful somersault is a balanced landing which in turn requires that the performer 
finish the somersault with a pre-specihed body attitude. In this chapter, I explore 
attitude control techniques for producing stable landing configurations for the tucked 
somersault of a 3D biped running robot. 

The tucked somersault is distinguished by maintaining the tuck position during 
most of the maneuver. Humans 'tuck' by holding their knees close to the chest 
with the knee joints flexed to fold the lower legs under the body. An important 
dynamic feature of the tuck somersault is that in most humans it involves rotation 
about the major principal axis of inertia. Since a rigid body rotating about its major 
principal axis is directionally stable, we might expect the tuck somersault to be stable 
in the sense that the axis of rotation will tend to maintain its inertial orientation. 
This stability in turn simplifies the control of body attitude at landing in the tuck 
somersault. 

Directional stability of the spin axis during a tucked somersault means that a 
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Figure 3-1: Photograph of the 3D Biped used for experiments. The body is an 
aluminum frame on which are mounted hip actuators and computer electronics. Each 
ball and socket hip joint has three low friction hydraulic actuators that can position 
the leg forward and aft, side-to-side, and can rotate the leg along the axis of the leg. 
A hydraulic actuator within each leg changes its length, while an air spring makes 
the leg springy in the axial direction. Sensors measure the lengths of the legs, the 
positions of the hip actuators, pressure in the hip actuators and contact between 
the foot and the floor. Gyroscopes measure the inertial attitude of the body. An 
umbilical cable connects the machine to hydraulic, pneumatic, and electrical power 
supplies. Control computations are done my microprocessors, some located on board 
and some nearby in the laboratory. Communications cables connect all processors. 

somersault control strategy need deal primarily with avoiding over-rotation or under- 
rotation about the somersault axis. Rotation rate about the somersault axis can be 
controlled by changing rotational inertia. Since angular momentum must be conserved 
during flight, increasing inertia will slow down the somersault rate while decreasing 
inertia will increase somersault rate. If the performer knows the time of flight, then 
control of somersault rate provides a means of controlling the somersault angle at 
landing. 

In this chapter, I discuss somersault experiments with a 3D biped robot, (Figure 3- 
1). Figure 3-2 shows a sequence of photographs of the 3D Biped taken while the robot 
performed a successful somersault. The somersault axis of the robot is coincident 
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with the maximum principal axis ol inertia making the robot somersault dynamically 
similar to the tuck somersault in humans. To initiate the somersault the biped runs 
forward, jumps to attain a double stance phase, then thrusts with both legs while 
pitching its body forward. Once airborne, the robot shortens its legs (tucks) to 
accelerate the forward rotation and swings its legs to a predetermined position with 
respect to the body. During flight the robot uses a feedback algorithm that changes 
the leg length to produce a rotation rate that will yield the desired somersault angle 
at landing. To accommodate errors in the estimated time of landing, the robot 
moves its feet to track the desired landing configuration as the system approaches the 
ground. The robot does not use any active control of out-of-plane rotation during 
flight. Rather it uses a broad stance during takeoff to minimize tilt rotation at the 
beginning of flight and it uses a broad stance during landing to minimize the effect 
of landing tilt errors. The robot has successfully performed the somersault in the 
laboratory. On its best day, the robot regained balance on landing to continue stable 
running on seven out of ten somersault attempts. 

3.2 The Mechanics of the Somersault 

The most basic requirement of a somersault is that the performer neither over-rotate 
nor under-rotate the landing. Accurate control of the landing attitude allows careful 
placement of the foot relative to the center of mass of the robot which is a requirement 
for stable dynamic running [Raibert 84]. Considering only the somersault degree of 
freedom, the attitude requirement is expressed by equating the time of flight and the 
time to rotate through the desired change in somersault attitude, 



A$ _ z + \jz 2 + 2g(z - z td ) 

$7 " * (3 ] 



where 
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Figure 3-2: A sequence of photographs (arranged in left-right order starting upper 
left) taken during the execution of a somersault. The robot is running from left to 
right. Approximate relative time of each image: upper left-0.0 s, upper right-0.15 s, 
middle left-0.33 s, middle right-0.66 s, lower left-0.80 s, lower right-1.02 s. 
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A$ is the change in somersault angle required during flight 

$ is the somersault rate of the body at lift-off 

z is the height of the center of mass (c.o.m.) at lift-off 

z is the vertical velocity of the c.o.m. at lift-off 

z t d is estimated height of the c.o.m. at touchdown 

g is the acceleration of gravity 

Equation 3.1 relies on several simplifying assumptions: 1) the somersault dynamics 
are governed by the planar equation /$ = 0, implying that the external torques due to 
supply cables or wind resistance are negligible, and 2) only the rotation in somersault 
is significant and the somersault axis is a maximum principal axis of inertia so that 
tilt and twist angles will stay small if they start small thus allowing us to ignore them, 
and 3) the legs do not swing with respect to the body during the flight phase, so $ 
represents the angular rates of both the body and the legs. 

When the 3D Biped robot somersaults it rotates about its major principal axis. 
Figure 3-3 shows that the region of stable rotation about the somersault axis is large. 
We may then expect that as long as the somersault is initialized with the angular 
momentum vector close to the major principal axis then it will remain close to that 
axis. This in turn means that the tilt and twist angles of the robot at landing will 
be small and the somersault dynamics simplify as indicated above. In the remainder 
of this chapter we assume that this simplification is valid in computations involving 
rotational dynamics of the robot. 

3.3 Somersault Control Strategies 

The goal of the somersault control strategy is to produce a landing attitude that allows 
the robot to maintain balance. To regain balance on the landing, it is important that 
the robot achieve a desired horizontal displacement of the foot relative to the center 
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Figure 3-3: The somersault axis of the 3D Biped is coincident with the major principal 
axis of inertia. With the legs in the fully extended position, as shown here, the minor 
principal axis is the 'head-to-toe' axis. 
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of mass. In the plane of the somersault this horizontal displacement is given by 

d = /sin/3 (3.2) 

where / is the leg length and f3 is the leg inclination angle, the angle the leg makes 
with vertical in the plane of rotation (Figure 3-4). If the legs are held fixed relative to 
the body, then the desired landing attitude, or equivalently A$, can be found from d. 
If the initial conditions of flight, $ , i , and z , are established accurately so that for 
a desired z t d and A$ Equation 3.1 is satisfied, then the desired landing attitude will 
be achieved. More generally, one approach for generating somersaults is to establish 
accurately the initial conditions of flight to a state that is empirically determined, 
then to execute an open loop pattern of actuator signals to produce components of the 
desired behavior. Hodgins and Raibert used this approach in programming a planar 
biped robot to perform front somersaults with a 90% success rate. The success of 
such an approach depends upon how precisely one can reproduce the state of the 
robot and how sensitive the desired movement is to variations in the state. 

For a running robot with less regular and repeatable motion, such as the 3D 
Biped, precise initialization is more difficult. Therefore, reliable production of a 
desired landing attitude might be improved with an in-flight feedback strategy that 
modifies the performance of the somersault based on the state of the robot in flight. 

To the extent a system with non-zero angular momentum can change its inertia, 
it can also change its rotation rate. For a somersault, if the time until landing is 
known then control over rotation rate amounts to control over the landing attitude. 
The basis of the somersault control strategy is to change the robot somersault inertia 
to produce a rotation rate that will yield the desired body attitude at the time of 
landing. 

The robot can change its inertia by extending or retracting (untucking or tucking) 
its prismatic legs. In flight, if the angular rate and moment of inertia of the robot in 
one configuration are $i and I\ then with constant angular momentum the angular 



CHAPTER 3. ROBOT TUCKED SOMERSAULTS 



56 




Figure 3-4: In order to maintain balance at landing the leg inclination angle, /3, must 
be near the desired value. f3 is the angle the leg makes with the local vertical. It is a 
function of the the hip angle, a and the body attitude, $. 

rate in a configuration with inertia I 2 is $ 2 = (hfh)®i- Joint limits restrict the 
range of inertia and thus limit control over the rotation rate. For the 3D Biped 
robot, the inertia about the somersault axis ranges from 1.22 kg m 2 with legs retracted 
to 1.50kg m 2 with legs fully extended. In moving from a fully retracted to a fully 
extended position the robot can reduce its somersault rate to 81% of its initial value. 
We will refer to the regulation of somersault rate and landing attitude via leg length 
as the tuck servo. 

3.3.1 Control of Somersault Angle 

To implement the tuck servo with state feedback, we first pose the requirement of 
Equation 3.1 as a function of the state during flight rather than at lift-off. 



$tj- $ 



' z o + \fz 2 +2g(z - cos j3 u l td ) 



$ 



(i-io) 



(3.3) 



where 
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<&td is the desired somersault attitude at touch down, 

$ is the current somersault attitude, 

$ is the current somersault rate, 

f3 t d is the desired leg inclination angle at touch down, 

ltd is the leg length at touch down, 

t is the current time, 

t is the lift off time. 

The robot inertia and thus $ are functions of l t d } so we solve for the ltd that will make 
Equation 3.3 an equality. 

We nominally require the hip angle to be zero on landing. If this is the case then 
fttd = *&td an d the landing attitude and leg length determine the height of the hip. 
The center of mass is coincident with the hip, so z t d is given by 

z td = cos $ td ltd (3.4) 

We have also assumed that the tilt angle at touchdown is zero. 

Define a function, /, as the difference between the left and right hand sides of 
(3.3). 



, z + ^f z 2 +2g(z - cos $td ltd) $*<*-$ ,_ ., 

/= (t-to) t (3-5) 

9 $ 

If the desired landing attitude is to be achieved then / = 0. Otherwise, the pitch 

rate needs to be increased or decreased depending on the sign of /. The control 

algorithm we use changes the leg to a length that will make / = 0. The function, /, 

is a non-linear function of the leg length, so we use a Newton search [Strang 86] to 

recursively solve for the desired leg length. The recursion uses the first order, Taylor 

series expansion of / = 0, 

= /(Z) + ^AZ (3.6) 
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This yields the following equation for A/ 

A ' = wki (3 ' 7) 

Ideally, once the desired leg length is achieved, / remains zero for the rest of the 
somersault and the desired leg length becomes the leg length at touch down, / — > l t d- 
During each control cycle the tuck servo executes the following process 

1. measure the state of the robot 

2. compute / and df/dl t d 

3. compute Al t d 

4. estimate the next value of / using l t d = ltd + Al t d- 

5. if / ~ go to (6), else go to (2) using new values of / and l t d 

6. servo the leg length to l t d- 

To perform step 2 above we need to solve for df/dl t d- Taking the derivative of 
(3.5) with respect to l t d results in 

df/dhd = . - cos ^ + ^LZ± S/dhd (3.8) 

IP + 2g(z - cos $ td l td ) $ 2 



To obtain d<&/dl t d we use the fact that the angular momentum is a constant so that 

• h 

$ = (J + 2m,r*) (3 ' 9) 

where 

h angular momentum 

r distance from the lower leg c.o.m. to the robot c.o.m. 

mi lower leg mass 
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I robot pitch inertia about c.o.m. when r = 0. 

The total inertia of the robot as a function of leg length is (7 + 2mir 2 ) . The 2mir 2 

term represents that inertia due to the distance of the lower legs from the robot c.o.m. 

This is the component of inertia that we control as the legs change length. From (3.9) 

we get 

d$/dl td = - ' dr/dlu (3.10) 

(1 + mir 1 ) 1 

The robot center of mass moves very little when the legs are extended or retracted 
so we assume dr/dl td = 1. Substituting for h from (3.9) results in 

d$dl td =-——± 3.11 

(1 + mir 1 ) 

Now by substituting (3.11), (3.8), and (3.5) into (3.7), we can compute the required 
Al td . 

The computation of the desired leg length depends upon knowledge of robot pa- 
rameters such as inertia and leg mass. However, since the process is repeated each 
control cycle based upon the sensed state of the robot, sensitivity to precise knowl- 
edge of these parameters is reduced. In exchange for this robustness to uncertainty 
we give up the ability to pre-specify both the landing attitude and the leg length at 
touchdown. 

3.3.2 Accommodating Landing Time Errors 

A limitation of the tuck servo strategy is its dependence upon accurate knowledge 
of the time until landing. With somersault rates on the order of 500 deg/s, small 
errors in the predicted landing time can result in intolerable landing attitude errors. 
Because we have no measurement of the vertical position or speed while airborne these 
quantities must be derived from estimates of the lift-off conditions. To accommodate 
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errors in these estimates, we have modified the landing strategy to allow the robot 
to land with the desired configuration anywhere within a 'window' of the predicted 
landing time. 

The leg inclination angle in (3.2) is the difference between the somersault angle 
and the hip angle, f3 = $ — a. By changing the hip angle the control algorithm 
makes the foot track the desired displacement for a short time just prior to or after 
the expected landing time. In this way, the robot maintains the desired landing 
configuration during a 'landing window' that is centered about the nominal landing 
attitude. In order to maximize the landing window the legs are moved 'back' relative 
to the body (a > 0) early in the flip. Prior to landing, the foot is swept forward in 
order to track the desired foot position until touchdown. 

3.3.3 Control of Tilt and Twist Angles 

Achieving a takeoff attitude with no tilt and twist is challenging with the 3D Biped. 
As the robot runs its body somersault angle stays close to zero but its body tilts and 
twists in phase with the stepping cycle as it runs. During a normal running cycle, 
the robot uses hip actuators during stance to apply torques to the body in order to 
control body attitude. However, we found that it was difficult to achieve small tilt 
and twist angles and rates at takeoff using only the hip servos. We found that the 
best method for keeping tilt and twist angles and angular rates small during takeoff 
was to use a wide double stance during the flip initiation in combination with the 
normal attitude control used during stance ([Raibert 84]). This wide double stance 
helps stabilize the robot's tipping motion during takeoff. An important component of 
this approach was to achieve simultaneous touch down of both feet at the beginning 
of the double stance phase. Similarly, a wide double stance during somersault landing 
and simultaneous touchdown of both feet were observed to minimize the effect of tilt 
and twist errors on the somersault recovery. 
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3.4 Experiments with 3D Biped Somersaults 

The 3D Biped is a two-legged robot that is tree to translate and rotate in 3D space. 
Each leg is mounted adjacent to the center of gravity of the body with a ball and 
socket hip joint. The hip joint allows leg rotation about the x, y } and z axes 
(±20°, ±60°, ±15° see Fig. 1). Hydraulic actuators control each of these degrees 
of freedom. The robot's telescoping legs contain a fourth hydraulic actuator that acts 
in series with a pneumatic compression spring. Mass properties of the 3D Biped are 
included in Table 3.1 Some of the kinetic energy of the machine is stored in com- 
pression of the air spring during each bounce, and returned to power the subsequent 
flight phase. Energy is added to the hopping oscillation by actively compressing the 
air spring with the hydraulic piston during stance. The 3D Biped maintains balance 
while running by performing three control tasks [Raibert 84]: 

1. during stance, the robot maintains body posture by applying hip torques be- 
tween the legs and the body, 

2. during stance, the robot adds energy to the air spring to maintain the hopping 
oscillation, and 

3. during flight, the robot positions the foot in anticipation of the next stance 
phase in order to control forward velocity. 

To execute a somersault, the 3D Biped modifies three steps in an otherwise normal 
running sequence. The robot performs a hurdle step during which it hops higher than 
normal as it prepares to land on both feet for the flip step. The flip step is initiated 
by thrusting with both legs while pitching the body forward. During the landing step 
the robot lands on both feet then resumes a normal running gait. The control actions 
used to execute the flip are summarized in Table 3.2. 

In laboratory experiments, the 3D Biped has successfully performed the forward 
somersault and regained a stable running cycle afterwards (Figure 3-2). 
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Table 3.1 
3D Biped Parameters 



Symbol 


Description 


Quantity 


m 


total mass 


31.44 (kg) 


mi 


lower leg mass 


0.652 (kg) 


Io 


somersault inertia 


1.02 (kgm 2 ) 


'rain 


min. lower leg radius 


0.404 (m) 


''rain 


min. leg length 


0.647 (m) 


''max 


max. leg length 


0.862 (m) 



Three sets of data from a successful somersault and a nearly successful somersault 
are shown in Figures 3-5, 3-6, and 3-7. Figure 3-5 shows data for the approach, flip, 
landing, and continuation of running for a successful somersault. The nominal desired 
landing attitude was set to 350°. The robot ran steadily until the hurdle step at which 
time it hopped higher than normal as it prepared to land on both feet. During the 
flip step the body is thrust upwards and accelerated in somersault. The desired leg 
inclination angle was set to —5.7° based on the forward speed in flight. The robot 
lands 0.080s earlier than anticipated with a somersault attitude of 325° and with a 
leg inclination angle of —11°. The actual leg inclination angle is much closer to its 
desired value than the corresponding values in somersault because of the feet-back 
position in flight. The foot positioning servo was not used in this somersault because 
the desired foot position was always slightly in back of the actual foot position, and 
the legs were already in a swept back configuration. Balance is regained on landing, 
but since the robot lands with a slight backward lean, forward speed is lost. Forward 
speed and posture are quickly restored during the following steps. 

Note the oscillation in tilt angle during flight. This oscillation occurs because the 
robot took off with a non-zero tilt rate. Since rotation about the somersault axis 
is passively stable, the oscillation does not grow. We found that in order to regain 
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Table 3.2 
Control Summary for the Somersault 



Step 


Action 


Approach 


Run forward (ss lm/s) with alternating gait 


Hurdle 


Pitch up slightly 




Hop with increased leg thrust 




Prepare to land simultaneously on both feet 


Flip 


Jump with maximum thrust 




Pitch body forward with maximum torque 




Shorten legs once airborne 




Servo hips to feet back position 




Engage tuck servo 




Prepare to land simultaneously on both feet 




Track desired foot position 


Landing 


Dissipate energy on landing 




Return somersault rate to zero, restore posture 




Adjust nominal leg length based on l t d 


Following 


Resume running with alternating gait 



CHAPTER 3. ROBOT TUCKED SOMERSAULTS 64 

balance after the somersault, the tilt angle on landing must be kept moderately small, 
1 1 < 15°. To do this we increased the stance width during the flip step, thereby 
providing a passively stable stance configuration in tilt during initiation. 

The data in Fig. 4 shows the somersault action on a larger scale to illustrate the 
function of the tuck servo. Between the beginning of the flip step and lift off, the 
magnitude of the somersault rate, hip angle rate, and leg length all increase. At lift 
off, the somersault rate of the body declines rapidly as the legs are accelerated to 
the rotation rate of the body. The conservation of momentum constraint produces 
the symmetry between the absolute angular rate of the body and the relative angular 
rate of the hips. During this time, the legs are tucked to the shortest possible length. 

The tuck servo is engaged as the hips reach the desired feet-back position and 
come to rest relative to the body. The robot has a somersault rate of 606 deg/s at 
the time the tuck servo is engaged. At this somersault rate, it is estimated that the 
robot will over-rotate by 84°. This error is illustrated by the third graph of Fig. 4 
which compares the estimated time until touchdown and the estimated time until 
the desired somersault attitude is achieved. The tuck servo extends the legs to the 
maximum possible length to slow down the somersault rotation to 463 deg/s at which 
point it is estimated that the robot will land at nearly the desired attitude. The robot 
maintains this configuration until landing. 

Data from another nearly successful somersault is shown in Fig. 5. In this somer- 
sault, the robot uses the foot positioning servo to keep from over-rotating. Once again 
the desired somersault landing attitude was 350°. The desired leg inclination angle 
was set to —5.9°. The robot lands 0.065s later than anticipated with a somersault 
attitude of 391° and with a leg inclination angle of —3.2°. As the robot detected 
that it was over-rotating it swept the feet forward quickly to track the desired foot 
position. This increased the somersault rate because of the conservation of angular 
momentum and contributed to the somersault attitude error on landing. A velocity 
measurement error during stance after landing led to the loss of forward speed of the 
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Table 3.3 
3D Biped Flip Attempts 



File 


Outcome 


Ptd 


$td 


$ 


i (est) 


i (meas) 


z %err 


speed /o 






deg 


deg 


deg/s 


m/s 


m/s 




m/s 


92.346.4 


success 


0.12 


-10.8 


547 


3.08 


3.04 


-1.2 


1.04 


92.346.5 


success 


-0.14 


-6.29 


529 


3.48 


3.33 


-4.3 


0.77 


92.346.6 


success 


0.33 


-8.47 


524 


3.41 


3.32 


-2.8 


0.89 


92.346.7 


success 


-0.31 


-13.7 


549 


3.53 


3.48 


-1.4 


1.10 


92.346.8 


6 steps 


-1.95 


-18.9 


553 


3.69 


3.70 


1.6 


1.02 


92.346.9 


success 


-0.39 


-12.7 


551 


3.45 


3.47 


0.7 


1.04 


92.346.10 


success 


-0.04 


-11.6 


545 


3.54 


3.48 


- 1.7 


0.96 


92.346.11 


success 


-0.28 


-15.3 


541 


3.15 


3.13 


- 0.5 


0.90 


92.346.12 


fall 


-3.42 


-1.04 


529 


3.28 


3.04 


- 7.4 


0.89 


92.346.13 


fall 


-1.87 


-17.7 


503 


3.33 


3.20 


-3.9 


0.88 



robot so that it was momentarily supported by safety ropes before resuming running 
and therefore not considered a complete success. 

Table 3.3 shows data compiled for ten somersault attempts performed by the 3D 
Biped in the laboratory. The robot successfully regained balanced running on seven 
of these attempts. 



3.5 Summary 

In this chapter, I discuss a strategy for robot somersaults that combines elements of 
feed forward control, feedback control, and passive dynamic stability. I also presented 
results from somersault experiments done in the laboratory on a 3D biped running 
robot. The somersault is initialized using pre-programmed patterns of action. In 
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Figure 3-5: Hip height, somersault, and tilt angles, and forward speed of the 3D-Biped 
during a somersault. Vertical lines indicate initiation of flip, lift-off and landing. Data 
hie B92.181.3 
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Figure 3-6: Somersault angle, absolute somersault rate (solid) and relative hip rate 
(dashed). The third graph shows the estimated time until touch down (solid) and 
estimated time until desired somersault attitude is achieved (dashed). The last graph 
is the measured leg length (solid) and desired leg length(dashed) for one leg. Vertical 
lines indicate initiation of flip, lift-off, tuck servo initiation, and landing. Data hie 
B92.181.3. 
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Figure 3-7: This somersault was over-rotating so the feet were swept forward at the 
end to track a desired foot position. Data shows absolute somersault rate (solid) and 
relative hip rate (dashed) and estimated time until touch down (solid) and estimated 
time until desired somersault attitude is achieved (dashed). Vertical lines indicate 
initiation of flip, lift-off, tuck servo initiation, and landing. Data hie B92.181.4. 
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flight, a feedback strategy changes robot inertia to control the landing attitude of 
the somersault. The robot actively positions its feet to maintain a desired landing 
configuration during an interval surrounding the predicted landing time. The passive 
tilt stability inherent in a wide double stance is used to reduce tilt angle and rate on 
takeoff. The passive stability of a rigid body rotating about its maximum principle 
axis of inertia accounts for moderate tilt angles on touch down given moderate tilt 
angles and rates on lift-off. 



Chapter 4 

Passively Stable Layout 
Somersaults 



4.1 Introduction 

The layout somersault is an airborne maneuver in which the performer rotates about 
a side-to-side axis while maintaining an erect body configuration. The layout somer- 
sault is often considered to be inherently unstable because it involves rotation about 
the middle principal axis of inertia, an unstable mode of rotation for a rigid body. 
I found that passive arm movement could neutrally stabilize the layout somersault. 
These passive movements are generated by dynamic forces that arise from body mo- 
tion and from compliant shoulders. Figure 4-1 shows a sketch of the type of three-body 
model used for simulation, analysis, and experimental tests in the laboratory. This 
simple three-body model of a human can exhibit passively stable layout somersaults 
if an appropriately tuned spring is used at the shoulder. 

Layout somersaults are considered stable if steady somersaulting is achieved with- 
out exhibiting any 'large' oscillations in the tilt or twist angles. For example, Fig- 
ure 4-2 shows simulation data of the simple human model for two cases: 

1. Rigid body (dashed lines): this model has very stiff shoulder springs and 
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Figure 4-1: Diagram of a dynamic model used to study layout somersaults. The head, 
torso, and legs comprise a single rigid body. The arms are connected to the body 
with pin joints that allow rotation about the 1 body axis. Joint torques are provided 
by torsional springs and (possibly) dampers. 
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dampers to approximate the rigid body case. As with the true rigid body, the 
tilt and twist angles periodically move through large excursions from zero as 
the body somersaults. 

2. Passively stable body (solid lines): this model has tuned shoulder springs 
without dampers that allow considerable passive arm movement. The tilt and 
twist angles of the body stay near zero. 

The rigid body tilt and twist angles shown in Figure 4-2 are shown in another 
form in the map of Figure 2-7 which is repeated here in Figure 4-3. The rigid body 
instability is evident by the trajectories that converge and diverge upon the somersault 
axis. The twisting oscillation of the simulation data corresponds to a trajectory on 
the map that is centered about the maximum principal axes (front-to-back). 

In this chapter, I present results of studies on the passive layout somersault using 
non-linear dynamic simulation and linear stability analysis. By deriving a linearized 
model of the dynamics of a three-body model, I show that passive layout stability 
depends upon both the arm angle and the shoulder spring constant. I show condi- 
tions under which passive layout stability can theoretically be achieved. Using the 
simplest possible model of the layout somersault. I explain the fundamental dynam- 
ics of passive stabilization. I also present summarized results of non-linear dynamic 
simulations that support the linear stability analysis. 

4.2 Simple Human Model 

Figure 4-1 shows a simplified human body model used to study the dynamic stability 
of somersaulting motion. From the non-linear equations of motion of this model I 
analytically derive linearized equations that govern the motion of the body relative to 
pure somersaulting rotation. The linear equations of motion decouple into two distinct 
subsystems which allow further model simplifications. When inertial parameters for 
a human are used in this model, a linear stability analysis indicates that stability of 
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Figure 4-2: Number of complete somersault revolutions, twist and tilt Euler angles 
describing body attitude, and left and right shoulder angles for two simulations. The 
solid lines indicate the passively stable case. The dashed lines indicate the rigid body 
case. The data for the passively stable case shows nearly eight stable somersaults 
about the middle principal axis. 
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Figure 4-3: A map of rigid body rotation for a human in the layout position. The 
inertially fixed angular momentum vector (black) paints trajectories onto the surface 
of the sphere as the sphere rotates. The axes of the sphere and the principal axes of the 
body remain parallel as the body rotates. Each trajectory on the sphere corresponds 
to a different rotational energy. The trajectories indicate the tilt and twist angles 
of the body as it rotates. This map does not include somersault angle and does not 
show the time dependence of the tilt and twist Euler angles. 
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the layout somersault can be achieved for a range of nominal arm angles and shoulder 
spring constants. 

The human body model has shoulder pin joints that allow the arms to be raised 
and lowered in the '2-3' plane of the body. The head, torso, and legs are modeled as 
a single rigid body. Torques at the shoulder joints are provided by torsional springs 
and dampers. Neglecting translation of the center of mass, the model has five degrees 
of freedom (d.o.f.): three for rotation of the body with respect to inertial space and 
two for relative movement of each arm. The orientation of the body with respect 
to an inertial coordinate frame is described using a 2-1-3 sequence of Euler angles 
(Figure 2-1.) To describe the attitude of the body relative to the inertial frame, 
a coordinate system initially parallel to the inertial reference frame is first rotated 
through the somersault angle about the inertially fixed '2' axis, then rotated through 
the tilt angle about the intermediate '1' body axis, and finally rotated through the 
twist angle about the body fixed '3' axis. 

The rotation of the right and left arms relative to the body are given respectively 
by the angles *y r and 7/. An arm angle is equal to zero when the hand is held next 
to the thigh. Positive rotation of the arm about the '1' axis is given by the right 
hand rule. The right and left arm angles are rewritten in terms of symmetric (7 S ) 
and asymmetric (^ a ) components as follows: 

Is = 1/2(7,-7,) 

7a = l/2( 7/ + 7r ) 

The parameters necessary to describe the three-body system in the non-linear 
equations of motion are given below. The body is assumed to be symmetric from left 
to right when the arms are held in a symmetric orientation. 
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Figure 4-4: Diagram of a five d.o.f. model with body axes and vectors labeled. 
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shoulder spring and damping constants 



Figure 4-4 shows vectors and body axes defined in the above list. 

In order to study the layout somersault as performed by humans, I will use data 
for a human performer, Carl Furrer, 1982 World Trampoline Champion [Yeadon 84]. 
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Table 4.1 
Inertia Data for Carl Furrer 



Parameter 


Value 


Units 


m b 


58.3 


kg 


mi 


3.743 


kg 


IW 


10.19 


kg m 2 


Ib 2 


9.93 


kg m 2 


Ib 3 


0.52 


kg m 2 


Ik 


0.128 


kg m 2 


Ik 


0.128 


kg m 2 


Ik 


0.0037 


kg m 2 


rb\ 


0.0 


m 


rb 2 


0.18 


m 


rb 3 


0.479 


m 


rl\ 


0.0 


m 


rl 2 


0.0 


m 


rk 


0.262 


m 



The data for Furrer 's body parameters are included in Table 4.1. 

4.2.1 Nonlinear Equations of Motion 

I derive the non-linear equations of motion 1 as a starting point for a linear analysis 
that will follow. I used Kane's method [Kane 85] to derive the equations. Rather than 
include the lengthly equations here I include the Mathematica code used to produce 
them in Appendix A.l. The equations are derived using a reference coordinate frame 
located at the center of mass of the three body system. Therefore, the rotational 
equations are decoupled from the translational equations. 



1 Non-linear equations used for dynamic simulation were derived using a commercially available 
package, SD-FAST. This version of the equations was not amenable to analytic linearization. 
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4.2.2 Linearized Equations of Motion 

The linearized equations will describe the motion of the system with respect to a 
reference frame that is steadily rotating about the 2 axis with rate $. Pure somersault 
rotation about a principal axis of inertia is an equilibrium solution for this system if 
a constant feed forward shoulder torque is applied to cancel the centrifugal forces due 
to steady rotation. This steady torque can be achieved by pre-tensioning the shoulder 
spring. The magnitude of this torque will be a function of the nominal symmetric 
arm angle, T s0} and the rotation rate, $. Any deviation of the arm from T s0 will 
result in additional shoulder torques from the springs and dampers. 
The linearized equations of motion for this model follow: 



Mx + $ (G + D) x + $ 2 (K + K') x = 



(4.1; 



where 



M 



G 



ran mis 

m 22 

m 13 m 33 

m 44 

m 55 






912 











912 





923 











~923 























945 











-945 






(4.2) 



(4.3) 
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(4.4) 



(4.5) 



(4.6) 



2k sh /& 

The state vector, x, is comprised of (in order): 6 - the tilt angle, ip - the twist 
angle, ^ a - the asymmetric deviation of the two arms from the nominal arm angle, 
<f> - the deviation of the somersault angle from the frame steadily rotating at rate $, 
and 7 S - the symmetric deviation of the two arms from the nominal arm angle. 

The expressions for individual components of Af, G, and K are listed in Appendix' A. 3. 
In Equation 4.1 the nominal rotation rate, $, scales the rate terms and $ 2 scales the 
spring terms. This is evidence of the fact that the coriolis forces involving somersault 
rate give rise to rate dependent terms and the centrifugal forces involving somersault 
rate give rise to spring terms. The terms in D and K' arise from the spring and 
damper model of the shoulder torques (77, r r ) for each arm. 



n = -k sh (-fi - r s0 ) - b sh ii 



(4.7) 
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T r = -k s h{^fr + T s0 ) - b sh i r (4.8) 

Non-dimensional states (z,i,£) are derived by scaling i, and x by $ or $ 2 as 
follows 

Z = X 

z = i/$ 
z = 5/$ 2 

Rewriting the equations of motion as functions of the non-dimensional states results 
in the following 

$ 2 (Mz +(G + D)z + (K + K') z) = (4.9) 

It is now clear that $ serves only as a time scale of Equation 4.9. We may consider 
$ = 1 for simplicity in studying stability as long as we are willing to work with 
scaled shoulder spring and damping terms. In the following stability analysis we use 
a non-dimensional shoulder spring defined as follows 



$ 2 L 



where I s h is the inertia of the arm about the shoulder. For Furrer's data I s h = 
0.3849% m 2 . 

The reason for describing arm angles in terms of symmetric and asymmetric com- 
ponents is clear upon examination of the matrices of Equations 4.2 - 4.6. This de- 
composition allows us to see that in general the linearized system decouples into two 
distinct subsystems; the upper 3-by-3 system that couples tilt and twist with asym- 
metric arm movement and the lower 2-by-2 system that couples rotation rate with 
symmetric arm movement. We consider these two subsystems separately in a linear 
stability analysis. It will turn out that the parameters which stabilize the composite 
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system are nearly, but not exactly, identical to those that stabilize the 3-by-3 system. 

4.3 Stability Analysis 

In this section, I present a linear stability analysis of the layout somersault. To start, 
I present the stability of the equivalent rigid body, that is the rigid body that would 
result if the arms were rigidly fixed to the body. As expected, if the body were rigid 
the layout somersault would be unstable. Then I proceed to an analysis of the three- 
body system. The goal is to determine the shoulder spring constants and arm angles 
that produce stable somersaults when the remaining body parameters are fixed. A 
root locus plot shows how the location of the linear system poles move as the shoulder 
spring constant is systematically varied. The root locus plot shows that under some 
conditions, all poles of the linear system are simultaneously on the imaginary axis, 
implying (neutral) stability. The results of the root locus are also presented in the 
form of stability diagrams that show under what values of shoulder spring and arm 
angle the passive stability can be realized. 

4.3.1 Rigid Body 

As a start, let us examine the characteristics of the equivalent rigid body for this 
system. I consider the equivalent rigid body to be the rigid body that has the same 
principal inertia as the whole three body system. If we consider the arms to be fixed 
rigidly to the body then we can solve for the equivalent rigid body inertias, and thus 
the inertia ratios, k\ and k 3} (see Section 2.3.2) as a function of the nominal arm 
angle, r s0 - The inertia ratios k\ and k 3 capture the inertia of the equivalent rigid 
body and therefore determine its stability. 

For a true rigid body rotating about its intermediate axis, the unstable mode of 
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motion will grow at a rate given by Equation 2.13 

z = z e s where s = (-k^) 11 ' 2 ^ 

Since $t is simply the nominal somersault angle, the change in somersault angle, $, 
required for the unstable mode to grow by a factor of N is given by Equation 2.14 



Using the data for Furrer (Table 4.1) in the three d.o.f. model with a symmetric 
arm angle of T s0 = 1.3 rad. the resulting rigid body inertia ratios are k\ = 0.65 and 
k 3 = —0.87. Using these numbers with N = 10 results in $ = 3.1 rad, or the unstable 
mode will grow by a factor of 10 in less than one half somersault. Figure 4-5 shows 
a plot of the inertia ratios of the five d.o.f. model for symmetric arm angles between 
0.0 and 3.0 radians. These inertia ratios are superposed upon curves of constant $. 

4.3.2 Multi-Body 

Some useful information regarding the stability of this system is available by a simple 
examination of the form of the matrices in Equation 4.1. If shoulder damping, b s h, is 
zero then D = and Equation 4.1 describes a conservative gyric system. As discussed 
in Section 2.3.1, the strongest stability result one may expect from a conservative gyric 
system is neutral stability, i.e. the system will be a perfect oscillator. A sufficient 
condition for stability is for the net stiffness matrix to be positive definite, K-\-K' > 0. 
This is not possible since K + K' is not even full rank for the five d.o.f. system. This 
happens because the system is unstable to perturbations in somersault rate as was 
the case for the rigid body system in Chapter 2. However, even the three-by-three 
subsystem of tilt, twist, and asymmetric arm movement does not have a positive 
definite stiffness matrix for the fixed body parameters used here. Therefore, if the 
system is to be stable it must be gyrically stabilized. 
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Figure 4-5: Rigid body stability diagram for Furrer's data. The dark lines show k\ 
and k 3 for < T s0 < 3.0. The light lines are lines of constant $, the change in 
somersault angle required for the unstable rigid body mode to grow by a factor of 
ten. The plot shows that the human form is relatively unstable in the rigid body 
sense. 
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Figure 4-6: Root locus plot for the 3-by-3 subsystem of tilt, twist, and asymmetric 
arm movement. b s h/<& = 0.0 and the nominal symmetric arm angle, T s0} is 1.3 rad. 
Light dots indicate root locations for K s h ~ —2.0. Dots get progressively darker as 
K s h increases to 4.0. 

Root Locus 



In this section I present plots of the roots of the characteristic equation of 4.9 as the 
shoulder spring constant and nominal symmetric arm angle are systematically varied. 
Figure 4-6 is a root locus for the 3-by-3 subsystem of tilt, twist, and asymmetric 
arm angle. This plot shows how the 3-by-3 system roots move as the shoulder spring 
stiffness is gradually increased when T s0 = 1.3 rad. Note the symmetry about the 
imaginary axis as predicted for a conservative-gyric system. 

When the nondimensional shoulder spring, K s h } is set to -2.0 the four symmet- 
rically located roots have both real and imaginary components. As the shoulder 
stiffness is increased the symmetric right and left half plane roots converge on the 
imaginary axis where they split. One set moves away from the origin along the imag- 
inary axis for as long as shoulder stiffness increases. The other pair converge at the 
origin where they again split one traveling along the positive real axis the other along 
the negative real axis. At very high shoulder stiffness these two real roots approach 
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the real (stable and unstable) rigid body roots. At high shoulder stiffness the high 
frequency imaginary pair correspond to an arm-body oscillation. 

There are two poles of this system at ±z that do not move with changing shoulder 
spring constant. These are the roots of the stroboscopic mode of rotation. This mode 
is a rigid body mode of motion that does not include arm movement. Examination 
of the eigenvectors associated with this mode reveals that the attitude and velocity 
vectors are parallel. This indicates that this mode involves pure spin about the middle 
principal axis, the equilibrium condition, when that axis is perturbed slightly from the 
original orientation. This rigid body mode will exist for all values of the parameters. 

For some values of shoulder spring all of the roots of this system are located on 
the imaginary axis. These marginally stable conditions suggest that for these values 
of the parameters the physical system may be stable. Note that even though equal 
size steps in K s h were used to numerically evaluate the root locus, the loci do not 
move in even steps. In particular the loci 'jump' to the imaginary axis from the right 
and left half planes. This indicates that the edge of the stable region in parameter 
space is steep. 

Figure 4-7 shows a root locus for the same system with the addition of a viscous 
damper at the shoulder joint. Note the loss of symmetry about the imaginary axis 
and the lack of a set of spring values that stabilize the system. This result agrees 
with the spirit of Hughes proof that damping tends to destabilize gyrically stabilized 
systems. However, this model does not exactly fit the requirements of his proof since 
the damping matrix is not positive definite. Finally, Figure 4-8 shows how the root 
loci plots change as T s0 is varied between 0.0 and 2.8 radians. 

The root loci of this section suggest that neutral stability of the layout somer- 
sault is possible for certain combinations of nominal arm angle and shoulder spring 
constant. It appears that asymptotic stability is unlikely as shoulder damping tends 
to destabilize the system. In the next section we will examine more precisely the 
dependence of stability on arm angle and shoulder spring constant. 
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Figure 4-7: Root locus plot for the 3-by-3 subsystem of tilt, twist, and asymmetric 
arm movement. b s h/<& = 0.1 and the nominal symmetric arm angle, T s0} is 1.3 rad. 
Light dots indicate root locations for K s h ~ —2.0. Dots get progressively darker as 
K s h increases to 4.0. 

Stability Diagrams 

In this section, I present plots that show precisely how stability depends on the 
shoulder spring constant and nominal arm angle. 

Figure 4-9 shows the results of a search in the T S Q-K s h parameter space for condi- 
tions that stabilize the five d.o.f. linearized system. This collection of points suggests 
a region in the parameter space that stabilizes the system. The stable points in Fig- 
ure 4-9 are an intersection of the stable points of the two decoupled linear systems 
that comprise the five d.o.f. system. Figures 4-10 and 4-11, are the corresponding 
plots of stable points for the 3-by-3 subsystem of tilt-twist-asymmetric arm and the 
2-by-2 subsystem of somersault-symmetric arm. Comparison of Figures 4-9 and 4-10 
shows that stability of the five d.o.f. system is well represented by the stability of 
the simpler 3-by-3 system, i.e. few conditions that stabilize the 3-by-3 system fail to 
stabilize the 5-by-5 system. 
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Figure 4-8: Root loci for the 3-by-3 subsystem of tilt, twist, and asymmetric arm 
movement as the nominal symmetric arm angle, T s0} is varied between 0.0 (top left) 



and 2.8 rad (bottom center) , b 



sh 



0, -2.0 < K sh < 4.0. 
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Figure 4-9: Each darkened point in this plot corresponds to a specific choice of K s h 
and the nominal symmetric arm angle (T s0 ) that stabilize the full five d.o.f. linear 
system, b s h/<& = 0. 
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Figure 4-10: Each darkened point in this plot corresponds to a specific choice of 
K s h and the nominal symmetric arm angle (T s0 ) that stabilize the three d.o.f. linear 
system of tilt, twist, and asymmetric arm movement, b s h/<& = 0. 
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Figure 4-11: Each darkened point in this plot corresponds to a specific choice of K s h 
and the nominal symmetric arm angle (T s0 ) that stabilize the two d.o.f. linear system 
of somersault and symmetric arm movement, b s h/<& = 0. Higher values of K s h than 
are shown on this plot are suspected to also stabilize the system. 
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4.3.3 Non-Linear Stability via Dynamic Simulation 

The linear stability analysis (Figure 4-9) predicts that passive layout somersault sta- 
bility can be achieved for nominal symmetric arm angles between 0.8 and 2.0 radians 
(assuming only positive shoulder spring constants are allowed.) We would normally 
expect the linear analysis results to be valid only in some region of the linearizing 
condition. Stability to perturbations away from the linearizing condition will depend 
upon the size of the perturbations. Since in the case of passive layout somersaults 
the linear analysis can only predict neutral stability we can not conclude that the 
non-linear system is stable even for arbitrarily small perturbations. While non-linear 
dynamic simulation can not prove non-linear stability it can provide a fast and easy 
check on whether or not the non-linear system is likely to be stable. In this section, 
I use non-linear dynamic simulation as a check on the linear stability results. 

Figure 4-12 shows results summarizing a series of dynamic simulations of the five 
d.o.f. human model. Each three dimensional plot shows how many stable somersaults 
the five d.o.f. model exhibited in simulation as a function of the nominal symmetric 
shoulder angle, T s0} and the shoulder spring, K s h- Each point in the grid represents 
a separate simulation. A somersault was considered stable if it did not exhibit a 
twist angle of at least ±7r/2 rad. There was little ambiguity between stable and 
unstable layout somersaults. Either the twist angle remained small while the model 
somersaulted or the twist angle would grow significantly beyond one-quarter twist 
making it easy to distinguish stable and unstable motion. The plots show that the 
region in parameter space that the linear analysis predicted to be stable is also the 
region that produces the most stable layout somersaults in non-linear simulation. 

The height of the plots in Figure 4-12 has been limited to eight somersaults. 
Several simulations exceeded this limit. For selected cases, the non-linear dynamic 
simulations indicated that the maneuver may remain stable indefinitely as long as 
energy is conserved. Two simulations starting from an initial tilt angle of 0.01 radians 
and zero twist angle were stopped after running for 2000 simulated seconds. During 
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these simulations the performer produced over 2750 stable layout somersaults. The 
shoulder springs and nominal arm angles that produced these long running maneuvers 
were {K s h,T s o} = {1.2, 1.4}, {1.0, 1.6}. Other conditions may also have produced 
stable, long running maneuvers. I did not perform an exhaustive search. 

4.4 How Passive Layout Stability Works 

In the last section we showed that under certain conditions a five d.o.f. dynamic 
model of a human can exhibit passively stable layout somersaults. We used numer- 
ical searches to find combinations of shoulder springs and nominal arm angles that 
stabilized the layout somersault. This search process does not provide much insight 
into how or why passive stabilization works. In this section, I describe what I consider 
to be the fundamental dynamic processes that make passive layout stability possible. 
First, I look at the differences in state trajectories between the passively stable case 
and the unstable rigid body case. The differences in trajectories between the two 
cases make evident that arm motion can change the orientation of the principal axes 
in a way that stabilizes the layout somersault. We then look more closely at how 
this arm movement arises. To do this I focus on the simplest possible model that 
captures passive layout stability. I define the important non-dimensional parameters 
that govern this model and show how stability depends upon these parameters. Then 
I provide an intuitive explanation of the passive dynamic process that stabilizes the 
layout somersault. 

4.4.1 Stabilization Via Principal Axes Reorientation 

Batterman [Batterman 68] and later Frohlich [Frohlich 79] showed how an athlete can 
control the orientation of his principal axes during a rotating maneuver using arm 
movement. This control allows an athlete to initiate twist in a somersault without 
using any external forces or torques. Using symmetry, Yeadon [Yeadon 84] has argued 
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Figure 4-12: These plots show that the regions (in T s0 - K s h space) of stable som- 
ersaulting motion of a non-linear dynamic simulation are similar to those regions 
predicted to be stable by the linear analysis (Figure 4-9). The height of each grid 
point indicates the number of stable layout somersaults performed during a separate 
dynamic simulation of the five d.o.f. model. The axes of each grid indicate the sym- 
metric arm angle, 0.4 rad < T s0 < 2.2 rad } and the non-dimensional shoulder spring, 
0.0 < K s h < 1.8, used for each simulation. The different plots show that the number 
of stable somersaults performed is sensitive to the initial tilt and twist angles. 
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Figure 4-13: The orientation of the principal axes of inertia can be controlled with 
arm rotation. The tilt of the principal axes of the figure on the right is a result of 
body axis tilt due to airborne reorientation of the body plus principal axes tilt due 
to body asymmetry. 



that an athlete has an equal opportunity to remove twist from a twisting somersault 
using arm motion. In this section I will show how arm movement is used to control 
the orientation of the principal axes. Then I will show how principal axis orientation 
can be used to influence rotational maneuvers. 

The orientation of a human's principal axes relative to a set of body fixed axes 
depends upon the configuration of the arms and legs of the performer. Figure 4-13 
shows a human figure in two different body configurations. The tilt of the body on 
the right is the result of the airborne reorientation of the arms from the symmetric 
configuration of the figure on the left to the asymmetric configuration of the figure 
on the right. Also shown are inertially fixed axes and principal axes of the figures. 
The principal axes of the figure on the left are aligned with the inertial axes. Notice 
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that the principal axes of the figure on the right are tilted with respect to the inertial 
axes and also inclined relative to the body axes. This inclination of the principal 
axes from the body fixed axes is due to asymmetry of the body and it has the affect 
of changing the 'gain' on principal axis reorientation due to arm movement. To see 
how principal axis orientation can be used to control twist we will examine a plot of 
rigid-body rotational trajectories. 

Figure 4-14 shows the twist angles plotted as a function of tilt angle for several 
rigid body trajectories of the human model while in the layout configuration. This 
plot is essentially a close up view of the {ci, c 3 } axes of the map shown in Figure 4-3. 
The different trajectories correspond to different energy levels of rigid body rotational 
motion. While time is not explicitly shown on the plot, the data points are all 0.01 sec. 
apart. Thus dot spacing shows that twist rate increases with tilt magnitude. Arrows 
have been drawn on this plot to show the direction of motion. The instability of 
rotation about the middle principal axis is evident from the converging and diverging 
trajectories near the origin. These trajectories show that if body attitude deviates 
slightly from the origin then it will depart the origin in a direction that depends upon 
which 'side' of the equilibrium it is located. 

Tilting the principal axes using arm movement has the effect of moving the body 
from one energy trajectory to another. While an instantaneous change in tilt can 
not change the twist angle, it can change the twist rate by putting the body on a 
trajectory with a different twist rate. Stabilization of the layout somersault can be 
viewed as a process of selectively adding or subtracting twist over time using arm 
movement so that the trajectory remains in the vicinity of the origin of Figure 4-14. 
Figure 4-15 shows a trajectory of a non-linear dynamic simulation of the passive, five 
d.o.f. model of the human. This trajectory of tilt and twist angles is superposed upon 
the rigid body trajectories of the same model. The passively stable trajectory circles 
the origin in a counter-clockwise fashion periodically increasing then decreasing the 
twist and tilt angles. The counterclockwise rotation of the trajectory means that 



CHAPTER 4. PASSIVELY STABLE LAYO UT SOMERSA ULTS 



96 



* 


;,. 


' • . 




0.10 
0:05 




. • 


■ ' 


• 


1 


i * 


1- 


"\ 


.... 


| 


f *'" f V 


• i 


i 


0.3 


-0.2 


/T 


/ 


.... 


\ 


K \\ 


0.2 


0.3 
Twist (rad) 


• 


• ' 


. • * 


. • 


-0.05 
-0.10 


r . * 


• • . 


• . 


. 



Figure 4-14: Tilt angles plotted as a function of twist angle for several different energy 
levels of rigid body rotation of a human. Arrows indicate the direction of motion with 
increasing time. Dots are all 0.01 seconds apart. The arms of the rigid human model 
were abducted to a symmetric angle of 1.3 radians. 
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Figure 4-15: Tilt angles plotted as a function of twist angle for a non-linear dynamic 
simulation of the five d.o.f. passive model of the human. The simulated trajectory 
is superposed upon rigid body trajectories of the same model. This simulation was 
initialized with a tilt angle of 0.01 rad, zero twist angle and symmetric arm angle of 
1.3 rad. 



while the body twists 'with the flow' in the upper and lower halves of the plane it has 
to move 'upstream' with regard to the tilt angle. How is this accomplished? 

The system uses arm movement to accomplish this 'upstream' tilting action. In 
Figure 4-15 we see that by starting with an initial condition of positive tilt angle, zero 
twist angle, and arms symmetric the body begins to twist in the negative direction. 
From this initial condition the tilt of the principal axes apparently decreases despite 
the rigid body tendency to increase the tilt angle. This decrease in principal axis tilt is 
a result of the arms moving asymmetrically (left arm up, right arm down) in response 
to the negative twist rate. The arm tilting (rotation about the '1' axis) arises from 
the combination of somersaulting (rotation about the '2' axis) and twisting (rotation 
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Figure 4-16: A simple model that captures the essential dynamics for passive layout 
stability. The single arm is allowed to pivot about its center. The arm and body 
centers of mass are coincident. 

about the '3' axis.) This gyroscopic affect will tend to move the arm in a direction 
that reduces the tilt angle of the principal axes as long as the body is twisting. 

This oscillation is a passive behavior that emerges from the mechanical properties 
of the system interacting with the environmental forces due to rotation. The gyro- 
scopic forces which provide the inherent correction of this model depend upon the 
shape and position of the arm. To see how this passive behavior emerges we will look 
at a simplified model of the dynamics. 



4.4.2 The Simplest Model 

The simplest model that I can think of that captures the essential dynamics of passive 
layout stability is shown in Figure 4-16 This model has the feature that only asym- 
metric arm motion is allowed. Recall from Section 4.3.2 that symmetric arm motion 
could be ignored when looking for stability of the five d.o.f. model. Also, since the 



CHAPTER 4. PASSIVELY STABLE LAYO UT SOMERSA ULTS 



99 



center of mass of the arm is coincident with the center of mass of the body, the 
equations of motion will be significantly simplified as compared to the more complex 
human model. 

The linear equations of motion for this four d.o.f. model will show that the 
somersault degree of freedom again decouples from the remaining dynamics. Ignoring 
this degree of freedom and assuming no damping at the pin joint of the arm results 
in the following three d.o.f. linear model. 



Mx + $Gx + $ 2 (K + K') x = 



(4.10) 



M 



G 



K 



mn rai3 

m 22 

m 13 m 33 

512 

-Ql2 #23 

-g 23 

hi &13 

k 22 

&13 k 33 




K' = 

k sh i& 

where the state vector, x, is comprised of (in order) tilt, twist, and arm angle and 
k s h is the spring between the arm and the body. The model parameters are given in 
Table 4.2. 

This dynamic model still requires us to define the six inertia parameters of the 
body before we consider stability. Recall from Chapter 2 that the rotational dynamics 
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Table 4.2 
Simplest Model Parameters 



Parameter 


Definition 


mn 


Ih + Ih 


m 13 


Hi 


m 22 


Ib 3 + Il 3 


m 33 


Ih 


912 


(Ih + Ih) - (Ih + Ih) + (Ib 3 + 7/ 3 ) 


923 


-7/i + Ih - Ik 


hi 


(Ih + Ih) - (Ih + Ik) 


ki3 


Ih - Ik 


k 2 2 


-(Ih + Ih) + (Ih + Ih) 


k 3 3 


Ih - Ik 



of a rigid body depend only upon two non-dimensional inertia parameters. Following 
the spirit of this result, in the next section we reduce the model to non-dimensional 
form in order further simplify the analysis. 

Non-Dimensional Linear Equations of Motion 

I now present a non-dimensional version of the linear dynamic equations of motion of 
the model shown in Figure 4-16 for the case of steady somersaulting about the '2' axis. 
The purpose of this step is to discover the important non-dimensional parameters of 
this model. The states are made non-dimensional through scaling by the nominal 
somersault rate. Casting the equations in terms of non-dimensional parameters is 
accomplished by dividing each equation by the corresponding inertia term on the 
diagonal of the mass matrix. A few more simple algebraic manipulations will bring 
the system of equations into the following form. (In this form these equations do not 
reflect the symmetry that we expect of a conservative gyric system. However, they 
can easily be made to fit this canonical form.) 



$ 2 (MS + Gi + Kz) 







(4.1i; 
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The five independent non-dimensional parameters are defined as follows: 



h 

h 

J-rel 

kh 

Ksh 



(Ib 2 + Il 2 ) - (Ih + Ih] 
Ih + Ih 

(ih + ih)-(ih + ih: 



ih + ih 



ih 



Ih + Ih 
Ih - Ih 


Ih 

ksh 



J/l$ 2 



Equation 4.11 depends upon only five non-dimensional parameters. The param- 
eters k\ and k 3 are the principal inertia ratios for an equivalent rigid body system. 
They describe the inertia properties of the system that would result if the arm were 
fixed in place. Notice that if the arm degree of freedom were eliminated from Equa- 
tion 4.11 then the remaining equations are identical to the rigid body equations (2.11). 
I re i defines how big the inertia of the arm is relative to the inertia of the rigid body. It 
is important in determining how much the body moves in response to arm movement. 
The next term, kh is an inertia ratio analogous to k\ involving only arm parameters. 
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Table 4.3 
Furrer's Non-dimensional Parameters 



Parameter 


Value 


h 

J-rel 

kh 


0.65 
-0.87 
0.06 
-0.99 



It is important in determining the tendency of the arm to tilt in response to twist 
rate. 

We use data from the human model (Table 4.1) to derive sample values for these 
dimensionless parameters. These parameter values are shown in Table 4.3. These 
parameters are computed assuming the arms are held straight out to the side of the 
body. To compute the arm inertia of the model in Figure 4-16 I let Il\ = Il 3 = 
2.0 X 0.384:9kg m 2 , the sum of the arm inertias about their shoulder axes in the given 
configuration. 

With these numbers in mind, we can inspect the individual terms of the matrices 
of Equation 4.11 to help us develop an understanding of the relevant dynamic forces. 
Without arm movement the system reduces to the rigid body equations of motion 
which are unstable for the set of parameters in Table 4.3. For this reason I will 
focus on the arm related terms only. First consider the 3-2 term of G. This term 
is nearly equal to 2.0. If the twist rate of the model is negative then this term will 
contribute to a positive acceleration of the arm angle which is precisely the corrective 
tendency that we observed in the plot of Figure 4-15. Additional arm accelerations 
come from the centrifugal forces associated with somersault rotation. These forces 
are represented by the 3-1 and 3-3 terms in the stiffness matrix, K. If K s h = then 
the 3-3 term of the stiffness matrix is negative. This reflects the fact that if the 
arm is deviated from its equilibrium position then it will experience centrifugal forces 
from somersault rotation that tend to drive the arm further away from zero. This 
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is because the centrifugal forces on either side of the arm are opposite one another 
due to asymmetry and the forces increase with deflection of the arm. Without a 
shoulder spring this term could overwhelm the helpful gyroscopic accelerations from 
the 3-2 term in G to destabilize the system. However, the shoulder spring offsets 
this destabilizing force. In fact, a rule of thumb for choosing the shoulder spring 
constant to stabilize the passive system is to choose the spring that exactly cancels 
the destabilizing term due to centrifugal forces. 

Performing a root locus search for stabilizing spring constants of this simple 3 
d.o.f. system results in the following range of stabilizing shoulder springs: 

0.98 < K sh < 1.08 

which approximately cancels kl\ = — 1. 

This rule for picking stabilizing shoulder springs generalizes to the more compli- 
cated model of the human. The most important difference in the human model is that 
the nominal symmetric angle of the arms can range from 0.0 radians (hands next to 
legs) to 7r radians (hands overhead). The non-dimensional shoulder spring constants 
that stabilize the human model are plotted as a function of nominal symmetric arm 
angle in Figure 4-17 (repeated from Figure 4-10). Also shown in this plot is the value 
of K s h that would exactly cancel the spring term due to centrifugal forces on the arm. 
Again, the region of stabilizing spring constants seems to be defined by the negative 
of this centrifugal force term. Finally, it is interesting to note that for a nominal 
symmetric arm angle of 7r/2 the non-dimensional stabilizing shoulder spring constant 
is approximately K s h = 1.0. 

Stability Diagrams for the Simplest Model 

The stability analyses performed thus far have all assumed that we have known values 
of the body inertia parameters. The small number of non-dimensional parameters of 
Equation 4.11 make it possible to study how linear stability depends upon these 
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Figure 4-17: Plot of the stabilizing K s h for a the 3-by-3 subsystem of the human 
model. The solid curve shows the value of K s h that will exactly cancel the spring 
term due to centrifugal forces on the arm. 
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inertia parameters as well. In this section I present two stability diagrams that show 
how linear stability depends upon the parameters £4, k 3} 7 re ;, and K s h- The stability 
diagrams consist of a sequence of plots of the {£4, £; 3 } parameter space. The darkened 
regions of these plots represent values of £4 and k 3 that are stable for the specified 
values of I re \ and K s h- 

Figure 4-18 shows that the value of K s h that stabilizes the largest region of the 
unstable fourth quadrant of the {k\ } k 3 } parameter space is K s h ~ 1.1. Figure 4-19 
shows that larger the value of I re \ the easier it is to stabilize the bottom two quadrants 
of the {£4, £; 3 } parameter space. Large I re \ corresponds to relatively larger arms. 

The parameter kl\ is a shape parameter for the arm that describes the tendency 
of the arm to tilt in response to twisting motion. Figures 4-9, 4-10, and 4-11 indicate 
that the effectiveness of the arm to stabilize the layout somersault decreases as the 
arms are moved away from T s0 = tt/2 (hands held straight out to the side.) I believe 
this change in effectiveness is due the change in the parameter kl\. To check this, I 
allowed the shape of the arm in Figure 4-16 to reflect the symmetric arm angle r s0 - 
Then I searched for the K s h that would stabilize the simple three d.o.f. model with 
fixed body parameters representative of those for a human. The stable T s0 - K s h 
configurations are shown in Figure 4-20. The shape of this plot is suggestive of the 
shape in the stability diagrams for the five d.o.f. model. An important difference 
between the diagrams for the simple three d.o.f. model and the five d.o.f. model is 
that the center of mass of the arm for the three d.o.f. model is coincident with the 
pin joint. I think this may lead to the extra stable conditions in Figure 4-20. 

4.5 Summary 

In this chapter I presented results from a study of the stability of several simple 
models of the layout somersault. A rigid body analysis of the layout somersault 
predicts that the maneuver is inherently unstable. A linear analysis of a passive three- 
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Figure 4-18: These plots show that 1.0 < K s h < 1.1 produces the largest region of 
stable conhgurations for rotation about the middle principal axis for the "simplest 
model". Each dark point represents a stable configuration in k\ - k 3 space for rotation 
about a principal axis. kl\ = —1, I re \ = 0.1 Each plot is for a different value of K s h- 
K sh = 0.0,0.5,0.8,0.9,1.0,1.1,1.2,1.5 and 2.0 for plots positioned from left to right 
and down the page. 
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Figure 4-19: These plots show that stabilization gets easier with larger arms. Each 
dark point represents a stable configuration in k\ - k 3 space for rotation about a 
principal axis. kl\ = —l,K s h = 1.0. Each plot is for a different value of I re \. I re \ = 
0.01,0.02,0.05,0.1,0.2,0.3,0.4,0.5 and 0.6 for plots positioned from left to right and 
down the page. 
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Figure 4-20: This plot shows that the stability of the "simplest model" can reflect 
the dependence on arm position that we saw in the five d.o.f. model (Figure 4-9.) 
Each dark point represents a stable configuration of the simple three d.o.f. model. 
The shape of the arm of the "simplest model" is changed with the value of T s0 to 
resemble the nominal symmetric arm angles of the five d.o.f. model. 
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body human model ol the layout somersault predicts that neutral stability can be 
achieved for a range of nominal symmetric arm angles and shoulder spring constants. 
The model is considered passive because the joint and body motion is the result 
of the passive dynamic interaction of the bodies and the environmental forces. No 
active control or sensing is used in the model. Non-linear dynamic simulation of the 
passive three-body model shows that the most stable configurations conform to the 
linear analysis results. The number of stable somersaults that the three body model 
performed in dynamic simulation before exhibiting a twist instability is sensitive to 
initial conditions. However, for selected cases, non-linear dynamic simulations of the 
layout somersault indicated that the maneuver may remain stable indefinitely as long 
as energy is conserved. 

In order to provide a more intuitive understanding of how passive stabilization of 
the layout somersault works, we considered the impact of principal axis orientation on 
layout stability. I showed how the natural tendency of the arms to tilt in response to 
twisting movement of the body provides a built-in correction to the twist instability. 
The arm tilt forces the principal axes of the system to move in a direction that 
compensates for tilt and twist errors. This built-in correction eliminates the divergent 
tendency of the system as long as the compliance of the shoulders cancels the unstable 
centrifugal forces on the arms. This effect is clarified by studying a very simple two 
body model of the layout somersault. This simple model also allows us to derive five 
non-dimensional parameters that are important in studying passive layout stability. 
I present stability diagrams that show how stability of the simplest model depends 
upon these parameters. 



Chapter 5 



Layout Somersault Experiments 



5.1 Introduction 

To find out whether or not passively stable layout somersaults are physically possi- 
ble, we built a somersaulting 'doll' that can be tested in the laboratory. The doll has 
springy shoulders that allow arm movement in the frontal plane. The non-dimensional 
inertia parameters of the doll show that it is dynamically similar to the five d.o.f. 
model of the human performer studied in Chapter 4. The goals of the experiments 
were 1) to determine if passive layout stability was physically possible, and, 2) em- 
pirically determine the shoulder spring constants that best stabilized the motion. 

To initialize the somersaults, we built a mechanical launching device that throws 
the doll into the air with angular rotation about the somersault axis. The launcher 
helped us to achieve consistency in the experimental conditions and helped to mini- 
mize the human influence of the experiment. The operation of the launcher and doll 
is similar in appearance to a gymnast swinging around a horizontal bar then releasing 
the bar to perform a multiple somersaulting dismount. 

In this chapter I describe the human-like doll and launching device used during 
experiments. Using models from the previous chapter, I perform an analytic stabil- 
ity analysis of the doll to determine under what conditions the doll is theoretically 
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stable in somersault rotation. Then, I describe the somersault experiments in which 
we repetitively launched the doll with different shoulder springs to determine which 
springs best stabilized the doll. During experiments we found that the doll could con- 
sistently perform at least three and one half stable layout somersaults during flight. 
The doll may have performed more layout somersaults if we had been able to observe 
longer flight times. The laboratory ceiling height limited our opportunity to observe 
more somersaults. Also, we found good agreement between theory and experiment in 
finding the shoulder springs that best stabilized the maneuver. 

5.2 Experimental Apparatus 

In this section I describe the human-like doll used during experiments. I compare three 
non-dimensional parameters of the doll and the human model studied in Chapter 4 
to show that the two models are dynamically similar. I also describe the launching 
device used to initialize the somersaults during experiments. 

5.2.1 Human-Like Doll 

I built a mechanical version of the five d.o.f. human model studied in Chapter 4. This 
doll is a fifteen inch tall, scaled version of the human performer whose anthropometric 
data were used in the stability analysis of Chapter 4. A photograph of this doll is 
shown in Figure 5-1. The body is made of wood. The shoulder joint axes are 3/16 
in. diameter steel shafts that are fixed to the body via aluminum brackets screwed to 
the front and back of the doll body. The arms are made of 1/4" threaded steel rod. 
The arms are attached to the shoulder shaft via an aluminum block fitted with ball 
bearings that allow free shoulder rotation in the frontal plane. Three shoulder pulleys 
with 1/2 in., 3/8 in., and 5/16 in. radii are mounted to the aluminum block to be 
concentric with the shoulder shaft. Linear springs attached between the body and 
the shoulder pulleys act as torsional springs about the shoulder shaft. The different 
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Figure 5-1: Photograph of the mechanical doll used for experiments with layout 
somersaults. 

radii pulleys allow a single linear spring to be used as three different torsional springs. 
Two sets of springs were used in opposition to each shoulder. The shoulder springs 
could be pre-tensioned so that the arm had any desired equilibrium angle. The doll 
had a 1/4 in. diameter vertical shaft attached to its head. This shaft is used as a 
handle for the mechanical launcher to grasp. One 30g threaded weight was added to 
each arm. Moving the position of this weight on each arm allowed us to change the 
inertia of the arm without changing the total mass of the body. 

The physical parameters of the five d.o.f. doll without added arm weights and 
with the arm weight located at a radius of 0.043m from the shoulder center of rotation 
are given in Table 5.2.1. The arm inertias are given in an axis system located at the 
center of mass of the arm and parallel to the body axis system when the arm is held 
straight down at the side, $ s0 = 0.0. Data for this chapter are all for the doll with 
arm weight at a radius of 0.043 m from the shoulder. 
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Table 5.1 
Doll Body Parameters 



Parameter 


No weight 




30g weight 
r = 0.043 m 




m b 


0.466 


kg 


0.466 


kg 


mi 


0.0275 


kg 


0.0575 


kg 


IW 


0.00411 


kg m 2 


0.00411 


kg m 2 


Ib 2 


0.00391 


kg m 2 


0.00391 


kg m 2 


Ib 3 


0.000354 


kg m 2 


0.000354 


kg m 2 


Ik 


0.0000559 


kg m 2 


0.0000789 


kg m 2 


Ih 


0.0000559 


kg m 2 


0.0000789 


kg m 2 


Ih 


0.000001 


kg m 2 


0.000001 


kg m 2 


rb\ 


0.0 


m 


0.0 


m 


rb 2 


0.0413 


m 


0.0413 


m 


rb 3 


0.0664 


m 


0.0664 


m 


rl\ 


0.0 


m 


0.0 


m 


rl 2 


0.0 


m 


0.0 


m 


rl 3 


0.0832 


m 


0.0572 


m 



Figure 5.2.1 shows a plot of the rigid body inertia ratios of the doll and of the 
human performer for symmetric arm angles between 0.0 and 3.0 radians. These are 
superposed upon curves that show how quickly the unstable somersault mode grows 
for a rigid body. Comparison of the inertia ratios shows that the doll has nearly 
the same rigid body inertia ratios as the human model. Figure 5-3 shows a plot of 
the non-dimensional parameter 7 i for both the doll and the human model for arm 
angles between 0.0 and 3.0 radians. Comparing these curves shows that the doll 
has relatively more massive arms than the human which makes stabilization slightly 



easier. 



5.2.2 Launching Device 

We built a launching device so that we could initialize the somersaults with consis- 
tency and with minimum direct human influence. A photograph of the launcher is 
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Figure 5-2: Rigid body inertia ratios of the doll (dotted line) and the human model 
(solid line) for arm angles between 0.0 and 3.0 rad. 
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Figure 5-3: Plot of I re \ as a function of nominal symmetric arm angle, T s0} for the 
doll (dotted line) and the human model. 

shown in Figure 5-4. The launcher functioned by accelerating the doll in rotation 
around a horizontal bar then releasing the doll at a fixed angle around the bar. The 
launcher consisted of three components: 1) a horizontal bar, 2) a bearing supported 
toggle clamp that was free to rotate around the bar, and 3) a torsion spring used to 
accelerate the clamp relative to the bar. 

To initialize the somersaults, the steel rod fixed to the head of the doll was clasped 
in the toggle clamp. The doll and toggle clamp would then be wound around the hor- 
izontal bar to tension the launcher spring. When the doll and toggle clamp were 
released they would accelerate around the horizontal bar. When the toggle clamp 
lever hit a trigger bar, the toggle clamp popped open, releasing the doll into a somer- 
saulting, parabolic trajectory. The position of the trigger bar around the horizontal 
bar could be changed so that the toggle clamp released the doll at an angle that pro- 
duced the desired flight trajectory. The toggle clamp had adjustments that allowed 
the doll to be launched with initial tilt or twist angles, however, for all experiments 
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Figure 5-4: Photograph of the mechanical launcher used to initialize the doll somer- 
saults in experiments. 
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Figure 5-5: The points in this plot indicate the values ol K s h that theoretically will 
stabilize layout somersaults in the doll. The cross-hairs indicate the spring constant 
that best stabilized the doll during experiments. 

we set the initial tilt and twist angles to be approximately zero. 



5.3 Theoretical Predictions 



Figure 5-5 shows the theoretically determined stabilizing values of K s h as a function of 
symmetric arm angle for the doll. Also shown in this figure is the spring constant that 
best stabilized the doll with a nominal arm angle of T s0 ss 1.57 during experiments. 
Note that this plot of the stabilizing K s h is very similar to those for the human in 
Figure 4-9. The most notable difference is the larger range of stabilizing K s h for the 
doll. This is due to the larger I re \ of the doll. 
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5.4 Description of Experiments 

The experiments consisted of a sequence of launches of the doll to determine which 
shoulder springs stabilized the doll in somersault rotation. For each different shoulder 
spring used on the doll, we launched the doll twenty times using approximately the 
same launch spring pre-tension. The toggle clamp released the doll from the bar at 
a somersault angle of approximately -| rad. We video taped each experiment. Using 
the video tape for analysis, for each launch we recorded the number of somersaults 
performed before the doll exhibited a twist angle of approximately -| (one quarter 
twist). If during a dismount the doll never exhibited this large twist angle then the 
dismount was considered to be completely stable. For comparison, we also performed 
twenty launches with the arms of the doll fixed rigidly in place. 

Recall from Chapter 4 that the nominal symmetric arm angle does affect the 
value of the stabilizing spring constant. During experiments we attempted to keep 
the nominal symmetric arm angle approximately equal to 1.6 radians by initializing 
the somersault with the appropriate arm angle and shoulder spring pretension. Large 
centrifugal forces on the arms during acceleration around the bar necessitated that we 
devise a method for fixing the arm position until the moment of release from the bar. 
In order to initialize the somersault with the desired arm angle we propped up the 
arms with thin steel rods held between the arms and small indentures near the feet of 
the doll. The rods were held in place by the shoulder spring pretension. As the doll 
rotated about the bar, centrifugal forces pushed the arms towards the feet helping to 
hold the steel rods and arms in place. Once free of the bar, the doll began rotating 
about its center of mass which tended to lift the arms off the steel rods allowing the 
rods to fall clear of the doll. The nominal arm angle was difficult to set precisely as 
it depended not only on initial position of the arms but also on the pretension of the 
springs necessary to hold the arms at this angle during rotation. 

Figure 5-6 shows the beginning of two separate experiments with the doll. 
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Figure 5-6: The photo on the left show the rigid doll with arms clamped in place. 
This doll exhibits the twist instability. The photo on the right shows the doll with 
flexible arms performing a stable layout somersault. 
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5.5 Experimental Results 

There was little ambiguity between stable and unstable dismounts. Either the twist 
angle remained small while the doll somersaulted or it would grow significantly beyond 
one-quarter twist making it easy to distinguish stable and unstable cases. The rotation 
rate available from our launcher and laboratory ceiling height of 4.3 m limited the 
doll to approximately three and one half complete somersaults from release to landing 
in the hands of our human 'catcher'. Figure 5-7 shows the average number of stable 
somersaults performed as a function of the six different shoulder spring constants 
tested. The average is computed over twenty sequential launches for a single shoulder 
spring. The standard deviation for each condition is shown with error bars. The 
average somersault rate was 16.6 rad/sec, and I s h = 3.021 • 10~ 4 iVm 2 . 

In addition to recording the number of stable somersaults exhibited during each 
launch we recorded whether or not the doll still appeared stable at the end of the 
maneuver (did not exhibit a quarter twist). If the doll had not exhibited a quarter 
twist by the end of the maneuver then it was considered to be a completely stable 
dismount. The percentage of completely stable dismounts as a function of K s h is 
shown in Figure 5-8. 

Of the six conditions tested, (five springs plus rigidly fixed arms) one shoulder 
spring value clearly outperformed the rest. The shoulder spring, K s h = 0.82, per- 
formed 18/20 dismounts without exhibiting the twist instability. This high percent- 
age of completely stable dismounts led to the very small variance in Figure 5-7. If we 
could have observed more somersaults, the doll may have exhibited a variance more 
in accordance with the remaining data. The variance in number of stable somersaults 
appears to increase with softer shoulder springs. While the doll became more erratic 
under these conditions, it never-the-less occasionally performed a completely stable 
dismount. At the other extreme, the doll with fixed arms reliably exhibited the twist 
instability at about one complete somersault from release. 

The experimental results contained in Figures 5-7 and 5-8 are in agreement with 
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Figure 5-7: Plot of the average number of stable somersaults performed during twenty 
launches of the doll as a function of the shoulder spring constant, K s h- The standard 
deviation is indicated by error bars. The data for the rigidly fixed arms is indicated 
by the infinite value of K s h- The nominal arm angle was approximately 1.6 radians. 
The best value of K s h = 0.82 falls withing the range of stability predicted by theory 
as indicated in Figure 5-5. 
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Figure 5-8: Plot of the percentage (of twenty) of completely stable dismounts as a 
function of the non-dimensional shoulder spring constant, K s h- The two test condi- 
tions on the right side of the plot exhibited zero completely stable somersaults. The 
best value of K s h = 0.82 falls withing the range of stability predicted by theory as 
indicated in Figure 5-5. 
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the theoretical results for stable layout somersaults. Figure 5-5 shows that for a 
T s o ~ 1.57, the stabilizing shoulder spring should have values in the range 0.7 < 
Ksh < 1-4. The best experimental value for the doll, K s h ~ 0.8 falls in this range 
with performance falling off on either side. Furthermore, this best value also falls 
in the range of stabilizing shoulder springs for the human model data as shown in 
Figure 4-9. 

5.6 Summary 

In this chapter I presented the results of layout somersault experiments on a mechani- 
cal human-like doll. Comparison of the non-dimensional model parameters of the doll 
and of a human indicate that the doll is dynamically similar to the human although 
slightly easier to stabilize due to relatively more massive arms. The experiments 
demonstrate that the doll can consistently perform at least three and one half sta- 
ble layout somersaults. The consistency of the dismounts depends strongly upon the 
value of the shoulder spring. The best shoulder spring exhibited 18/20 completely sta- 
ble dismounts, somersaulting dismounts without any evident twist instability. While 
softer shoulder springs make the dismount more erratic, they never-the-less allow the 
doll to occasionally perform a completely stable dismount and regularly perform two 
stable layout somersaults. On the other hand, stiff springs make the doll more reliable 
but less stable. The best value of the non-dimensional shoulder spring constant was 
K s h = 0.82 which is within the region of stability predicted by the linear stability 
theory. 



Chapter 6 



Twisting Somersaults 



6.1 Introduction 

The twisting somersault is a maneuver in which the performer simultaneously ro- 
tates about the somersault and twist axes of the body. Multiple twisting, multiple 
somersaulting maneuvers are among the most exciting and complex aerial maneuvers 
performed by gymnasts and other athletes. Unlike the pure somersault, the twisting 
somersault must include non-linear rotational coupling between the different body 
axes. One effect of this non-linearity is that the effect of the performer's control 
actions will change during a maneuver making cause and effect relationships more 
complex than in the pure somersault. Navigation and feedback control of twisting 
maneuvers are challenging tasks. Does the accurate, reliable performance of this 
maneuver necessitate a feedback control strategy? 

In this chapter I briefly discuss the open loop control of twisting somersaults. Open 
loop control means the performer's control actions are simply replayed from memory. 
We would like to know if twisting somersaults could exhibit passive dynamic stability 
when performed with an open loop control strategy. We discuss one test we performed 
to look for evidence of passive dynamic stability in the twisting somersault. 

I present results from the non-linear dynamic simulation of a I 1/2 twisting front 



124 



CHAPTER 6. TWISTING SOMERSAULTS 125 

somersault. This maneuver requires that the performer execute a sequence ol limb 
motions during flight. We found that when a simulated performer used a prescribed 
set of motions for executing this maneuver, the landing attitude of the performer was 
sensitive to initial conditions leading to poor landing attitudes. However, when the 
performer's control movements were compliant, the reliability of the landing attitude 
was significantly improved. There appears to be an optimal choice for the performer's 
compliance that leads to the most reliable landings. 

I also discuss the control of twisting somersaults by the 3D Biped robot. We 
programmed a simulated 3D Biped robot to perform a front somersault with half 
twist. We found that in order to make the robot maneuver look like a front somersault 
with twist as performed by a human, we had to add weight to the robot to make its 
moments of inertia more like those of a human. We also had to use stronger actuators 
than available for the physical robot. Our experiments to make the real 3D Biped 
robot perform the maneuver in the laboratory were unsuccessful. The physical robot 
actuators had insufficient actuator power to perform the maneuver. 

6.2 The Tilt of Twisting Somersaults 

Frohlich [Frohlich 80] described two techniques for performing twisting somersaults. 
In the torque twist the athlete derives rotation about the somersault and twist axes 
from external forces as he or she leaves the ground, diving board or other apparatus. 
In the torque-free twist with angular momentum the athlete initiates twist from 
an airborne somersault with an asymmetric movement of the limbs (Figure 1-1). The 
net effect of either technique is to tilt the principal axes of the body relative to the 
angular momentum vector. A map of rotation, such as that in Figure 2-7, shows that 
even a small amount of tilt of the principal axes from the layout somersault position 
will result in twisting. (This is what makes layout somersaults challenging.) The 
greater the tilt angle, the greater will be the twist rate. 
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6.3 One and One Half Twisting Front Somersault 

Yeadon [Yeadon 84] discusses a "torque-free" twist technique based on the hula move- 
ment for initiating twist from a piked front somersault. The hula movement involves 
a swiveling of the hips not unlike that required to swing a hula hoop about the hips. 
A quarter cycle of hula movement performed during a pike front somersault will tilt 
the performer's principal axes relative to the angular momentum vector. This effect 
is increased if the arms are held in an abducted position during the movement. After 
the hula movement the performer extends from the pike and the arms are adducted 
to decrease the inertia about the twist axis. The extension from the hula movement 
should happen between the 1/4 and 3/4 twist positions to maximize the net tilt of 
the body. However, between these limits, Yeadon claims the timing of the extension 
is not critical to the final tilt angle. 

6.3.1 The Nominal Case 

I created a simulation of a human performing a one and one half twisting front somer- 
sault. The simulated athlete used the technique described in the previous paragraph 
to perform the maneuver in a weightless environment (Figure 6-1). The human model 
has thirteen joint degrees of freedom, including three in each shoulder, one at each 
elbow, two in each hip, and one in the torso allowing the upper body to twist relative 
to the lower body. The maneuver was initialized from an upright piked somersault 
position with the arms held straight out from the side of the body. The maneuver 
finished in a layout body position with the arms held straight out from the sides. 

The control movements used to produce this maneuver were hand programmed 
so that the model looked natural during the movement and finished with the desired 
attitude and body configuration. The control movements consisted of a sequence of 
desired positions of the joints written as functions of time. The purpose of designing 
this maneuver was to empirically find a sequence of body configurations and joint 
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Figure 6-1: Images arranged in right-to-left, top-to-bottom order from a dynamic 
simulation of a 1 1/2 twisting front somersault. The maneuver was initialized from 
somersault rotation in the piked position. The control movements for this maneuver 
are a hand programmed sequence of joint angles written as functions of time. 

torques that would produce the desired maneuver under nominal conditions. The 
nominal somersault rate at the beginning of the maneuver was 12.6 rad/sec and the 
maximum twist rate during the maneuver is 40.2 rad/sec. 

6.3.2 Off-Nominal Performance, Prescribed Motion 

If the simulated performer uses the prescribed set of control movements from a differ- 
ent set of initial conditions than the nominal case then the trajectory of the maneuver 
will change. How sensitive the maneuver is to variation in the initial conditions is im- 
portant to the reliability of the maneuver. We tested the sensitivity of this maneuver 
to initial conditions by running a series of simulations, each starting from a different 
set of initial tilt (0o) and twist (Wo) angles of the body. We evaluated the reliability 
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of the maneuvers by comparing the landing attitude of the off-nominal simulations 
to that of the nominal case. (The landing attitude was considered to be the attitude 
at a fixed time after the start of the maneuver.) The error, e in landing attitude was 
computed as follows: 

njoints 

e 2 = l/4(($ d - $) 2 + (Q d - 0) 2 + (M> d - M>f + 1/13 £ (qd, - qi f) 

8 = 1 

where $, 0, and *& are the body attitude at landing qi refers to the i th joint position 
and the subscript d refers to the desired value. This equation for the error emphasizes 
the body attitude over body configuration. An error of 1.0 is large; it could mean the 
twist angle or somersault angle was off by 4.0 radians or about 270°. 

The results of a series of simulations that varied the initial tilt and twist attitude 
of the body over a range, — 0.1 rad < O < 0.1 rad } —O.lrad < ^ < 0.1 rad are 
shown in Figure 6-2. 

The nominal maneuver, O = 0, \&o = 0, corresponds to the center grid point 
of this figure. The height of the surface there is zero. Away from the nominal the 
landing attitude error increases except for a narrow valley of initial tilt and twist 
attitudes along which the landing error remains small. It appears that a prescribed 
motion strategy would not produce reliable 11/2 twisting somersaults. Is there a 
simple open loop strategy that can improve this performance? 

6.3.3 Off-Nominal Performance, Compliant Motion 

It seems unlikely that people could accurately reproduce prescribed motions in a 
dynamic movement like the 11/2 twisting somersault. People use springy muscles 
and tendons to position their limbs. It seems likely that the change in environmental 
forces that would accompany a change in maneuver trajectory would mean that the 
limb movements change even if the athlete tried to execute the exact same motions. 
This idea is the basis for the following set of simulation experiments. 
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Figure 6-2: This figure shows the accuracy of a 1 1/2 twisting somersault in the 
presence of off-nominal initial conditions. The simulated performer used a prescribed 
set of joint angles to perform the maneuver. The two grid axes correspond to initial 
tilt and twist angles of the body, —0.1 rad < O < 0.1 rad } —0.1 rad < ^o < 0.1 rad. 
The height of the grid indicates the accuracy of the landing attitude of the simulation 
starting from the corresponding set of initial conditions. 
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We ran a series of simulations of the 11/2 twisting front somersault that used the 
joint torques from the nominal maneuver as feed forward commands. We also used 
the prescribed joint motion from the nominal maneuver as the commanded positions 
for a set of compliant, position plus derivative (pd) servos at the joints. In this 
open loop strategy, the net joint torques would be a combination of the feed forward 
torques plus pd servo torques. Since the pd servos use a prescribed set of positions 
as desired values, they act like passive springs and dampers at the joints. If the 
maneuver started from the nominal initial conditions then the joints would follow the 
prescribed trajectories. Then the torques from the pd servos would remain zero and 
the nominal maneuver would be reproduced. 

Our goal in this part of the experiment was to find a set of pd servo gains (spring 
constants and damping coefficients) that produced the most reliable performances. 
Rather than search over the gain parameters of pd servos at thirteen joints, I chose to 
search over one parameter. Therefore, I compute the gains of all joint servos according 
to a single parameter, the body clamped natural frequency, u n . The intent behind 
this choice is that all body joints have a similar compliance or natural frequency of 
operation. We will then look for a body natural frequency that produces reliable 1 
1/2 twisting somersaults. 

To compute the servo gains at a given joint as a function of io n} I assume that the 
body inboard (towards the torso) from the joint is inertially fixed and all out-board 
joints are immobilized. This way the model simplifies to a single d.o.f. joint between 
the limb in question and ground. The equations reduce to a simple second order 
system of equations as follows 

I x + c (i — Xd) + k (x — Xd) = 

where I is the apparent inertia of the limb at the joint, c is the damping coefficient, 
and k is the spring constant of the joint. This simple system can be rewritten in the 
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canonical form 

x + 2(u n (x - x d ) + to 2 n (x - x d ) = 

therefore c/I = 2C,io n and k/I = u\. We choose ( = 0.7 to achieve a well damped 
system response [Ogata]. Then choosing u n allows us to compute the stiffness and 
damping constants in a consistent manner. 

We systematically varied the value of u n between 150 rad/s and 20 rad/s. At each 
of these values of the body natural frequency we performed a series of simulations 
starting from initial conditions just as in the case of prescribed motion at the joints. 
I plotted the results in the form of the 3D plot shown in Figure 6-2 for each frequency 
and subjectively evaluated them. There was a clear choice for the best natural fre- 
quency. Figure 6-3 shows the simulation results for the case of u n = 30 rad/s. This 
plot shows that the reliability of the 11/2 twisting somersault performed with the 
open loop strategy shows marked improvement over the prescribed motion case. Now, 
nearly half the set of initial conditions results in small final attitude errors. Further- 
more, the landing attitude errors were worse for either smaller or larger values of 
u n . These basic results held true for variations in the angular momentum, h, of the 
maneuver as well (0.9 h < h < 1.1 h). 

Figures 6-4 and 6-5 show data from two simulations using the nominal value of 
angular momentum and with O = —0.1, and ^ = —0.1. The data of figure 6-4 
shows that when prescribed joint motion is used for an off-nominal maneuver, the 
body attitude error becomes large at the end of the maneuver. In contrast, the data 
of figure 6-5 shows that while the joint angles incur some error during the maneuver 
the body attitude is close to the desired value at the end of the maneuver. 

Since the timing of this maneuver will scale with somersault rate it is instructive 
to show how the body clamped natural frequency compares to the nominal somersault 
rate, ^p = ^|^ = 2.38. It appeared in simulations that the compliance that allowed 
the arms to open prematurely, thus increasing the body inertia about the twist axis, 
was important in off-nominal simulations. For this reason, one might consider scaling 
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Figure 6-3: This figure shows the accuracy of a 1 1/2 twisting somersault in the pres- 
ence of off-nominal conditions. The simulated performer used feed forward torques 
plus passive pd servos at the joints to perform the maneuver. Pd servo gains were 
chosen according to u n = 30rad/s. The pd servos used the prescribed joint angles 
from the nominal maneuver as desired values. The two grid axes correspond to initial 
tilt and twist angles of the body, —0.1 rad < O < 0.1 rad } —0.1 rad < ^ < 0.1 rad. 
The height of the grid indicates the accuracy of the landing attitude of the simulation 
starting from the corresponding set of initial conditions. 
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Figure 6-4: Simulation data from all/2 twisting somersault that used prescribed 
control motions for initial conditions O = —0.1, and ^ = —0.1 are shown with solid 
lines. The body attitude and joint angles of the nominal maneuver are shown with 
dashed lines. While the prescribed joint positions are accurate for the off-nominal 
case (bottom three graphs), the body landing attitude error is large. (The traces stop 
at the landing time.) 
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Figure 6-5: Simulation data from a 1 1/2 twisting somersault with u n = 30, ©o = 
—0.1, and ^o = —0.1 are shown with solid lines. The desired values of the body 
attitude and joint angles are shown with dashed lines. Joint angles incur significant 
tracking errors but the landing attitude closely follows that of the nominal maneuver. 
(The traces stop at the landing time.) 
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the body clamped natural frequency by the maximum twist rate of the nominal 
maneuver, -^- = -^- = 0.746. 

' * ml 40.2 

The body clamped natural frequency was chosen to be a simple parameter that 
described the natural frequency of the whole body. However, since the computation 
of this parameter assumed that parts of the body were inertially fixed, the actual 
eigenfrequencies of the system will differ from this body clamped frequency. The 
actual eigenfrequencies of the body for the initial configuration (pike) and the twisting 
configuration (wrap) of the Rudi, are shown in Table 6.3.3. 

6.4 Twisting Somersault of the 3D Biped 

We programmed a simulated 3D Biped robot to use the "torque twist" method to 
initiate a twisting somersault. To use this technique the simulated 3D Biped robot 
accelerated its torso in somersault and twist during the stance phase just prior to 
take-off. These two components of momentum should result in an angular momentum 
vector that is tilted with respect to the principal axes of the robot. 

Figure 3-3 shows the map of rotational motion of the 3D Biped. When compared 
to the map for the layout somersault in the human (Figure 2-7) we see that the robot 
requires significantly more tilt of the principal axes (relative to the angular momentum 
vector) to achieve a twisting somersault. This is partly due to the fact that the robot 
somersault axis is the major principal axis but also due to the significantly different 
shape of the regions on the two maps. The qualitative difference in shape of the maps 
is due to the difference in rigid body inertia ratios of the robot and of the human. 
To put it simply the human is long and skinny and the robot is short and fat. This 
difference makes the twist harder to achieve (more tilt required), and it means the 
maneuver will not look much like that of a human. 

In its current configuration, the twisting somersault mode of the 3D Biped would 
force the robot to lay on its side some time during flight. This body orientation is not 
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what we identify with a twisting somersault in a human. The human's long axis stays 
closer to the vertical during a twisting somersault. In addition, this inertial orientation 
of the robot is undesirable from a practical standpoint because this orientation is 
coincident with gimbal lock in the gyros used to measure body attitude. During 
laboratory experiments, the physical robot achieved this horizontal position which 
resulted in damage to the gyroscopes. To correct this situation we changed the robot 
inertia, and thus its rotational modes, to look more like that of a human. 

To change the robot inertia we added weight to increase the major and interme- 
diate principal inertias without increasing the minor principal inertia. We did this 
by adding weight along the '3' axis of the robot. Figure 6-6 shows the new map of 
rotation for the 3D Biped with a 4.0 kg weight added 0.8 m above the hips of the 
robot. This map looks much more like that of a human now and the twisting somer- 
sault should bear this resemblance as well. Figure 6-7 shows a sequence of computer 
graphic images of a simulated 3D Biped with added weight (and strong actuators) 
performing a front somersault with one half twist. Notice that the simulated robot 
lands the twisting somersault facing the opposite direction it started from. Data 
from the simulation is included in Figure 6-8. The simulated running robot regained 
balance on landing to continue stable dynamic running. 

We tried this maneuver with the physical 3D Biped robot in the laboratory. The 
added weight and inertia were too large for the robot to achieve sufficient flight time 
or angular momentum to produce the maneuver. 

6.5 Summary 

In this chapter I presented simulation results of a 1 1/2 twisting front somersault 
performed by a model gymnast with thirteen joint degrees of freedom. The maneuver 
is initiated from a piked front somersault. The twisting maneuver results from a 
sequence of movements of the limbs and torso made during flight. We found that 
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Figure 6-6: Spherical coordinate map of the tilt and twist Euler angle trajectories for 
the 3D Biped with 4 kg of weight added at a distance of 0.8 m above the hip. 
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Figure 6-7: Sequence of images of the simulated 3D Biped somersault with twist. 
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Figure 6-8: Hopping height, somersault, tilt and twist Euler angles of the 3D Biped 
during a simulated front somersault with 1/2 twist. The simulated robot passes 
through an Euler angle singularity at approximately 1.2 sec. causing discontinuties 
in the data. This singularity did not affect the dynamic simulation however as a 
different set of attitude parameters were used. 
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prescribed limb movements could produce reliable 11/2 twisting somersaults only lor 
a small set ol off-nominal initial conditions. In contrast, a control strategy that used 
feed forward joint torques in a tuned passive dynamic model ol the performer could 
produce reliable maneuvers for a much larger set ol off-nominal initial conditions. 

We do not yet know the significance ol this result. There are many degrees ol 
Ireedom ol this model. Too many to easily distinguish between important and unim- 
portant effects. Simplification ol this model may reveal salient features that can be 
analytically confirmed. These simulation results from the rudi are potentially inter- 
esting because it suggests that a passive dynamic approach to twisting somersaults 
may be capable of producing reliable maneuvers as was the case with the layout som- 
ersault. It is interesting that many ol the eigenlrequencies ol the most reliable, tuned 
compliant system for the rudi were in the neighborhood ol the maximum twist rate ol 
the maneuver. This appears to be analogous to the ideal choice ol arm-body oscilla- 
tion in the passive layout somersault. In this case the best shoulder spring resulted in 
an eigenlrequency ol the arm-body oscillation that was equal to the somersault rate. 

In this chapter, I also described simulation experiments with 3D Biped twisting 
somersaults. The distribution ol mass in the 3D Biped robot makes a twisting som- 
ersault particularly difficult. In order to enter a rotational mode that could produce 
the twist angles desired, the robot had to assume body attitudes unlike those asso- 
ciated with a human somersault with twist. This extreme body attitude could lead 
to mechanical difficulties with inertial instruments on the physical robot. We added 
weight to the robot so its natural rotational modes were more like those ol a human. 
A simulated 3D Biped robot was able to perform a 1/2 twisting Iront somersault. The 
simulated robot landed the maneuver and continued running stably afterwards. The 
added weight made twisting somersaults ol the physical 3D Biped robot impossible 
due at least in part to insufficient actuator power. 



Chapter 7 



Summary and Discussion 



Inspiration for this thesis comes from the remarkable performance of aerial maneu- 
vers by gymnasts, divers, trampolinists, and skiers. From a distance, the control of 
aerial maneuvers appears to be a very complex task. The shortage of control during 
flight and the non-linear rotational mechanics make this problem more than a little 
daunting. However, adversity breeds creativity and athletes and other people have 
been very creative in working with these difficulties to produce elegant solutions to 
movement control problems. It appears that we can learn something from human ca- 
pability in making dynamic flying or floating machines that can reorient themselves 
quickly, reliably, and accurately. 

This thesis is concerned with how people can incorporate known aerial movement 
techniques into a strategy that produces reliable performances. People can influence 
aerial maneuvers through the relative movement of their limbs and torso during flight. 
While a fixed sequence of control movements can be found to produce a desired 
maneuver, prescribed limb motion may provide little reliability when a maneuver is 
subject to the variety of conditions that would accompany any real performance. It 
is possible that people produce reliable maneuvers by sensing these variations and 
actively computing responses to compensate for them. However, the complexity of 
this approach warrants the search for a simpler strategy. A persistent goal of this work 
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has been to find control strategies that do not depend entirely on active feedback, 
that are simple, and that can work in laboratory machines. 

The contribution of this thesis is to point out via example that passive dynamic 
solutions to movement control are an alternative to certain problems previously con- 
sidered to require active feedback control. We focus on an open loop, passive dynamic 
strategy for producing gymnastic maneuvers. Open loop means that control effort 
(torques at joints) are simply replayed from memory during a maneuver (the motor 
tape model of biological control). Under identical conditions open loop control will 
produce identical maneuvers. Using open loop control under different initial condi- 
tions will produce different outcomes due to the inherent dynamics of the system. 
This may make open loop control by itself incapable of producing reliable maneuvers. 
However, it is possible that the passive dynamic behavior of the performer's body 
could automatically compensate for off-nominal conditions in a way that produces 
reliable maneuvers. The passive dynamic behavior of a system can be tuned through 
the selection of passive elements like springs (springy muscles) or by choosing nominal 
body configurations during the maneuver. Careful tuning of the passive dynamics of 
a system coupled with open loop control may provide a strategy for producing reliable 
gymnastic maneuvers. 

This strategy is simple in the sense that it requires no computation during the 
maneuver. Only memory playback is required. Compensatory movements are 'com- 
puted' by the physical system as part of its natural behavior. Incorporation of this 
strategy could simplify gymnastic performances by reducing the amount of active 
control and computation required by the athlete. 

To the extent possible, we have insisted upon laboratory experimentation to test 
our ideas. We feel that laboratory experimentation on physical machines forces one 
to discover the salient features of a problem that may be difficult to reveal otherwise. 
(Another reason is that laboratory experimentation is also a lot of fun!) However, 
while a laboratory experiment can be used to validate a strategy for machine control, 
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it can only suggest the viability of a strategy for human use. 

We studied three different gymnastic maneuvers in this thesis, the tucked somer- 
sault, the layout somersault, and the 11/2 twisting front somersault. 

7.1 Robot Tucked Somersaults 

The tucked somersault is a common maneuver in athletic events. The main require- 
ment of this maneuver is a balanced landing which in turn requires a precise body 
attitude. We programmed a 3D Biped robot to perform front somersaults in the 
laboratory. The control strategy relied implicitly upon the passive dynamic stability 
of the tucked somersault. We used an active feedback system to control somersault 
rotation rate and foot placement at landing. 

A linear analysis of rigid body rotation showed that the tilt and twist attitude of 
the somersaulting 3D Biped robot is passively stable. However somersault attitude 
is unstable. This well known result led us to use active control of the somersault 
degree of freedom in the 3D Biped maneuver. The robot controls somersault rotation 
rate by tucking or untucking its legs during flight to change the rotational inertia. 
During flight we depend upon a stable passive dynamic response of the tilt and twist 
dynamics to insure small tilt and twist angles at landing. To keep the initial tilt 
and twist angles small we use a wide double stance of the robot during take-off and 
insure that the feet touch down simultaneously. Another important element of active 
control of the 3D Biped somersault was the placement of the feet prior to landing. 
This element was necessary to compensate for errors in the estimated landing time. 
On its best day the robot has performed successful front somersaults on seven out of 
ten attempts in laboratory experiments. 

We do not know if active control of rotation rate in the 3D Biped robot is necessary 
to perform the front somersault. It is possible that open loop control plus active foot 
placement prior to landing could also be used to perform reliable front somersaults. 
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This was the case with the planar biped somersault work of Hodgins and Raibert. A 
strategy that we would like to experiment with in the future could provide a middle 
ground between the active rotation rate control used in this study and the (primarily) 
open loop approach used in the planar biped. This strategy would use tuned leg 
springs to passively influence somersault rotation rate. This technique could correct 
rotation rate errors by allowing the legs to extend at higher than normal rotation 
rates (due to higher centrifugal forces) and vice versa for slower rotation rates. This 
passive dynamic approach coupled with active foot placement may provide reliable 
front somersaults with less active feedback control. 

7.2 The Layout Somersault 

The layout somersault involves rotation about the middle principal axis of inertia, 
an unstable rotation for a rigid body. Biomechanics researchers have suggested that 
human athletes use active feedback control during flight to stabilize the layout somer- 
sault. We found that the layout somersault could passively be made neutrally stable. 

The layout somersault is stabilized by passive arm movement in the frontal plane. 
The inherent tendency of the arms to tilt in response to twisting movement of the body 
provides a built-in correction for the layout instability. The arm accomplishes this by 
changing the orientation of the principal axes with respect to the angular momentum 
vector. This built-in correction eliminates the divergent tendency of the system as 
long as a carefully selected shoulder spring is used to cancel unstable centrifugal forces 
on the arms. These results are confirmed with linear stability analysis, non-linear 
dynamic simulation and laboratory experiments with a somersaulting doll. During 
experiments the doll can consistently perform at least three and one half stable layout 
somersaults. 

Analysis of a simplified model of the layout somersault dynamics revealed several 
salient features of the movement. The simplified model revealed the simple rule 
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for picking stabilizing non-dimensional shoulder springs, K s h = 1.0. This spring 
constant insures that the natural frequency of the arm-body oscillation is matched 
to the rotation rate of the body. It also is the exact spring required to cancel the 
destabilizing centrifugal forces on the arms. The empirically determined, best value 
of K s h for the experimental doll was K s h = 0.82. This value is in agreement with the 
theoretical results. I think that it is significant that the stabilizing shoulder spring 
just cancels the spring-like action of centrifugal forces. This balance between opposing 
forces allows the arm to be responsive to the gyric forces which ultimately stabilize 
the sytem. 

The simple model also helped explain why arm orientation is important in stabi- 
lizing the layout somersault. The tendency of the arms to tilt in response to twist is 
greatest when the arms are held straight out to the side. This tendency is decreased 
as the arms are raised or lowered from that position. Another feature that the simple 
model revealed is that stabilization of the layout somersault is easier with bigger and 
bigger arms. This is not surprising for in the limit when the arms (held straight out 
to the side) are much larger than the body, the system resembles a rigid body rotating 
about its minimum principal axis, a well known gyrically stabilized configuration. 

We ignored the control of somersault angle in the layout somersault. Stabiliza- 
tion of the somersault angle will be important for practical layout somersaults as it 
was for the 3D Biped tucked somersaults. However, decoupling of the tilt, twist, 
and asymmetric arm movement from the somersault and symmetric arm movement 
suggests that a somersault control strategy and a tilt-twist control strategy could be 
developed separately then superposed in the complete system. 

7.3 Twisting Somersaults 

Twisting somersaults are among the most interesting and complex of aerial maneu- 
vers. They involve simultaneous rotation about the twist axis and the somersault 
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axis. Their dynamic complexity makes feedback control of these maneuvers chal- 
lenging. Is it possible that open loop control can be used to reliably produce these 
maneuvers? Perhaps it can if the performers body acts like a tuned passive dynamic 
system. Our results on the open loop, passive dynamic control of twisting somersaults 
are of a preliminary nature. However, they point to a potentially interesting subject 
for future work. These results suggest that passive dynamics may help the reliable 
performance of twisting maneuvers as was the case for the layout somersault. 

We used a thirteen joint human model to simulate the performance of a 1 1/2 
twisting front somersault. The twisting maneuver was initiated from a piked front 
somersault using asymmetric movement of the limbs. We found that using prescribed 
limb motion during the maneuver produced inaccurate landing attitudes when the 
initial body attitude was allowed to vary. On the other hand, an open loop, passive 
dynamic strategy was able to produce reliable landing attitude for a relatively large 
set of initial conditions. In this strategy, open loop torque commands from the nom- 
inal trajectory were combined with torques from passive springs and dampers at the 
joints. We tuned the passive compliance of the performer's joints to find the value of 
compliance that produced the most reliable maneuvers. We used a single parameter 
to characterize the compliance of all the body joints. There was a clear choice for 
the best system compliance. More or less compliance at the joints led to less reliable 
performance. Inspection of the eigenfrequencies of the model revealed that the most 
reliable system had natural frequencies in the vicinity of the maximum twist rate of 
the maneuver. As in the layout somersault, it may be that a balance between the 
centrifugal forces of twisting and the spring forces of the joints is required to make 
the maneuver reliable. 

While the reliability of the open loop, passive dynamic system was better than that 
of the prescribed motion system there is still room for improvement. There were sets of 
initial conditions that were 'close' to nominal for which landing attitude performance 
was seriously degraded. Could a different passive system improve reliability in this 
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region as well? Perhaps a simple active controller could complement the passive 
dynamics to produce an even more reliable system. Yeadon has discussed the need 
for good timing in certain maneuvers. Perhaps a simple active controller that changed 
the replay rate of a set of pre-recorded actions based on perceived errors could get 
the timing 'right' to further improve maneuver reliability. 

7.4 Do Humans Use Passive Dynamics 

We do not know if humans use a passive dynamic approach to stabilization of gym- 
nastic maneuvers. This thesis can not prove or disprove the human use of such a 
strategy. The results of this thesis can only suggest that such a strategy is a viable 
one for particular maneuvers, or parts of maneuvers. What experiments could we do 
to learn more about how people perform these maneuvers? 

It would be useful to measure the range of initial conditions from which athletes 
can reliably perform aerial maneuvers, and also measure their corresponding limb 
movements. If the variations in initial conditions and control movements were large 
then one may argue that dynamic compensation rather than a prescribed motion 
strategy was at work. One could look for correlations between different initial condi- 
tions and different control movements. Could a passive dynamic model explain any 
observed correlations? 

Tests of the compliance of human muscle could reveal if the effective spring con- 
stants of human limbs are near the values predicted by theory to produce passively 
stable maneuvers. This could support a passive dynamic theory but not prove it's use. 
Measurements of the electromyographic signal of the muscles may not help to decide 
what strategy people use. Even with accurate recordings of many performances of 
a maneuver it would be difficult to distinguish between a motor tape model and an 
active control model of muscle activation. 

A difficulty in making a distinction between an active control strategy and an 
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open loop, passive dynamic strategy is that both approaches could produce a simi- 
lar dynamic response. In fact, one may argue that if there is one clear technique to 
achieving a reliable maneuver then both strategies would have to use it. This possibil- 
ity is illustrated by the fact that the active controller designed by Yeadon to stabilize 
the layout somersault used the same control effect as the passive strategy presented 
in this thesis. Both techniques used arm tilt to control body twist. Yeadon's tech- 
nique depended upon sensing the twist rate of the body and using it to compute an 
appropriate arm tilt response. Whereas, in the passively stable case the arm tilt was 
produced from the inherent dynamics of the system. Passive dynamic control is a 
subset of the space of active controllers that require little if any on-line computation. 

7.5 A Passive Dynamic Theory of Control De- 
sign? 

Can passive dynamic control be formulated into a machine design and control the- 
ory? I did not use nor did I develop a consistent theory of passive dynamic design 
for controlling gymnastic maneuvers. However, beyond the educated guesses I used 
to find some answers, elements of linear and optimal control theory are central to the 
approach used to analyze the movements of this thesis. Linear analysis allows one to 
examine the behavior of a dynamic system in the neighborhood of a known solution. 
Biology provides us with many examples of movements that we know work. I have 
made explicit use of these known solutions to search for passive dynamic stability. 
Tucked and layout somersaults can be analyzed using linear methods. Linear analysis 
of the twisting somersault will be more difficult but may help reveal the important 
features of this maneuver. Searching for passive dynamic solutions to known move- 
ments could, I think, be developed into a theory. In fact, one may argue that a passive 
dynamic theory of control may consist of a restriction of the optimal control theory to 
passive dynamic elements. Extension of a passive dynamic theory to unknown move- 
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ment solutions will be more difficult, reflecting the complexity of non-linear systems. 
However, perhaps this is where the greatest pay-off from a passive dynamic theory 
could lie. 

Why is it so easy to tell animal movement from machine movement? The richness 
of human and animal behavior reflects the complexity of the dynamic systems, their 
bodies, they are compelled to use. But this richness does not just come from complex- 
ity, it comes also from coordination. Relatively speaking, machine movement is still 
in its infancy. It would be a grand goal of any designer to make a machine that could 
move like an animal. It is an unfortunate consequence of the complexity of non-linear 
systems that non-linear control design techniques frequently rely upon cancellation of 
the dynamics in order to reach a solution. This approach risks the design of forced, 
uncoordinated behavior. Perhaps we should search for design techniques that em- 
brace the inherent dynamics of a system rather than avoid them. Perhaps passive 
dynamic stability could act as a guiding principal in the design of graceful movement 
in machines. Perhaps it already serves this purpose in animal movement. 
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A.l Mathematica Code for Non-linear Equations 
of Motion 

Following is the Mathematica code for deriving the non-linear equations of motion of 
the five d.o.f. system. 

(* File for defining 3D equations of motion of a three link body with 
links connected separately to the central body by two pin joints 
allowing rotation about the x body axis. *) 

(* same as man5dof .m except arm angles defined in terms of symmetric 
and antisymmetric components phil = phia + phis, phir = phia - phis) 



(* inertia matrices *) 

ib = {{ibll,0,0},{0,ib22,0},{0,0,ib33}} 
il = {{illl,0,0},{0,il22,0},{0,0,il33}} 
ir = {{irll,0,0},{0,ir22,0},{0,0,ir33}} 



(* main body *) 
(* left limb *) 
(* right limb *) 



(* body to joint vector for left arm, in arm coord. *) 
bjl = {bjll, bjl2, bjl3}; 

(* body to joint vector for right arm, in arm coord. *) 
bjr = {bjrl, bjr2, bjr3}; 

(* inboard body to joint vector for left arm, in body coord. *) 
ibjl = {ibjll, ibjl2, ibjl3}; 

(* inboard body to joint vector for right arm, in body coord. *) 
ibjr = {ibjrl, ibjr2, ibjr3}; 
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(* transformation from body coordinates to limb left coordinates thru phi 
a pos. phi means limb left rotates positively relative to body *) 

cbtol = {{1, 0, 0}, 

{0, Cos[phia[t] + phis[t]], Sin[phia[t] + phis[t]]}, 
{0, -Sin[phia[t] + phis[t]], Cos[phia[t] + phis[t]]}} ; 

cltob = Transpose [cbtol] ; 

(* transformation from body coordinates to limb right coordinates thru phir 
a pos. phir means limb right rotates positively relative to body *) 

cbtor = {{1, 0, 0}, 

{0, Cos[phia[t] - phis[t]], Sin[phia[t] - phis[t]]}, 
{0, -Sin[phia[t] - phis[t]], Cos[phia[t] - phis[t]]}} ; 

crtob = Transpose [cbtor] ; 

(* inertia matrices of arms expressed in body coord *) 
ilb = cltob . il . cbtol; 
irb = crtob . ir . cbtor; 

(* vector from net eg to body eg, in body coord *) 

rhob = -(ml (ibjl - cltob. bjl) + mr (ibjr - crtob .bjr))/mt ; 

(* vector from net eg to left limb eg, in body coord *) 
rhol = rhob + ibjl - cltob. bjl; 

(* vector from net eg to right limb eg, in body coord *) 
rhor = rhob + ibjr - crtob. bjr; 
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(* angular velocity of body in body coord.*) 
wbi = {wbl[t], wb2[t], wb3[t]}; 

(* angular velocity of left limb w.r.t. body in body coord.*) 
wlb = {phia' [t] + phis' [t], 0, 0}; 

(* angular velocity of limb in body coord.*) 
wli = wbi + wlb; 

(* angular velocity of right limb w.r.t. body in body coord.*) 
wrb = {phia' [t] - phis' [t] , 0, 0}; 

(* angular velocity of right limb in body coord.*) 
wri = wbi + wrb; 

(* wbi tilde matrix for cross products *) 

wbitilde = {{0,-wbi [[3]] ,wbi [[2]] },{wbi [[3]] ,0, -wbi [[1]]}, {-wbi [[2]] ,wbi[[l]] ,0}}; 

(* wli tilde matrix for cross products *) 

wlitilde = {{0,-wli[[3]] ,wli[[2]]},{wli[[3]] ,0, -wli [[1]]}, {-wli [[2]] ,wli[[l]] ,0}}; 

(* wri tilde matrix for cross products *) 

writilde = {{0,-wri [[3]] ,wri [[2]] },{wri [[3]] ,0 ,-wri [ [1] ] } ,{-wri [ [2]] ,wri[[l]] ,0}}; 

(* velocity of body e.g. *) 

vbi = Expand [D [rhob,t] + wbitilde . rhob] ; 
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(* velocity of left limb e.g. *) 

(* all components of rhol in body coord, so cross with wbi *) 

vli = Expand [D [rhol ,t] + wbitilde . rhol]; 

(* velocity of right limb e.g. *) 

(* all components of rhor in body coord, so cross with wbi *) 

vri = Expand [D [rhor ,t] + wbitilde . rhor]; 



(* Matrix of Partial velocities : 

columns correspond to generalized speeds 

1 - wrt wbl[t] 

2 - wrt wb2[t] 

3 - wrt wb3[t] 

4 - wrt phia' [t] 

5 - wrt phis ' [t] 

rows correspond to body velocites 

1 - vbi 

2 - wbi 

3 - vli 

4 - wli 

5 - vri 

6 - wri 
*) 

pv = { 

{D [vbi, wbi [t]] ,D[vbi,wb2[t]] ,D[vbi,wb3[t]] ,D[vbi,phia' [t]] ,D [vbi, phis' [t]]}, 
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{D[wbi,wbl[t]] ,D[wbi,wb2[t]] ,D[wbi,wb3[t]] ,D[wbi,phia' [t]] ,D[wbi, phis' [t]]}, 
{D[vli,wbl[t]] ,D[vli,wb2[t]] ,D[vli,wb3[t]] ,D[vli,phia' [t]] ,D[vli, phis' [t]]}, 
{D[wli,wbl[t]] ,D[wli,wb2[t]] ,D[wli,wb3[t]] ,D[wli,phia' [t]] ,D[wli, phis' [t]]}, 
{D[vri,wbl[t]] ,D[vri,wb2[t]] ,D[vri,wb3[t]] ,D[vri,phia' [t]] ,D[vri, phis' [t]]}, 
{D[wri,wbl[t]] ,D[wri,wb2[t]] ,D[wri,wb3[t]] ,D[wri,phia' [t]] ,D[wri, phis' [t]]} 

>; 

(* inertia forces and inertia torques *) 

(* all vectors are expressed in body coordinates. See pp26 of notebook 
dated 6/4/93- for explanation of the following, ilb and irb are 
inertia of arms expressed in body coord therefore they are not 
constant in body coord, system. Therefore in taking total 
time derivative of ang. mom. need to take time derivative 
of ilb and wbi parts. Since everything in body coord, need 
to take cross product of ang. mom. with wbi. 

*) 



infb = mb wbitilde.vbi + mb D[vbi,t]; 

intb = D[ib,t] . wbi + ib . D[wbi,t] + wbitilde . ib . wbi; 

infl = ml wbitilde. vli + ml D[vli,t]; 

intl = D[ilb,t] . wli + ilb . D[wli,t] + wbitilde . ilb . wli; 

infr = mr wbitilde. vri + mr D[vri,t]; 

intr = D[irb,t] . wri + irb . D[wri,t] + wbitilde . irb . wri; 
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(* generalized inertia forces *) 



fil = pv[[l,l] 
pv[[3,l] 
pv[[5,l] 

fi2 = pv[[l,2] 
pv[[3,2] 
pv[[5,2] 

fi3 = pv[[l,3] 
pv[[3,3] 
pv[[5,3] 

fi4 = pv[[l,4] 
pv[[3,4] 
pv[[5,4] 

fi5 = pv[[l,5] 
pv[[3,5] 
pv[[5,5] 



infb + 


pv[[2,l]] 


infl + 


pv[[4,l]] 


infr + 


pv[[6,l]] 


infb + 


pv[[2,2]] 


infl + 


pv[[4,2]] 


infr + 


pv[[6,2]] 


infb + 


pv[[2,3]] 


infl + 


pv[[4,3]] 


infr + 


pv[[6,3]] 


infb + 


pv[[2,4]] 


infl + 


pv[[4,4]] 


infr + 


pv[[6,4]] 


infb + 


pv[[2,5]] 


infl + 


pv[[4,5]] 


infr + 


pv[[6,5]] 



. intb + 
. intl + 
. intr; 

.intb + 
.intl + 
. intr; 

.intb + 
.intl + 
. intr; 

.intb + 
.intl + 
. intr; 

.intb + 
.intl + 
. intr; 



(* External forces expressed in body coordinates *) 

(* 

positive acting on left limb, neg on body use wli,wbi 
phi_rest is the spring rest length, should be less than 
the nominal phil in order to make taul in steady rotation 
cancel the centrifugal force. 
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*) 



ftaul = {-ksh (phia[t] + phis[t]) - bsh(phia' [t] + phis' [t] ) ,0,0}; 
ftaur = {-ksh (phia[t] - phis[t]) - bsh(phia' [t] - phis ' [t] ) ,0,0} ; 



(*generalized external forces *) 

fel = -pv[ [2,1]] .ftaul + pv[ [4,1]] .ftaul 
fe2 = -pv[[2,2]] .ftaul + pv[ [4,2] ] .ftaul 
fe3 = -pv[ [2,3]] .ftaul + pv[ [4,3] ] .ftaul 
fe4 = -pv[[2,4]] .ftaul + pv[ [4,4]] .ftaul 
fe5 = -pv[[2,5]] .ftaul + pv[ [4,5] ] .ftaul 



pv [ [2 , 1] ] . ftaur + pv [ [6 , 1] ] . ftaur ; 
pv[[2,2]] .ftaur + pv[[6,2]] .ftaur; 
pv[[2,3]] .ftaur + pv[[6,3]] .ftaur; 
pv[[2,4]] .ftaur + pv[[6,4]] .ftaur; 
pv[[2,5]] .ftaur + pv[[6,5]] .ftaur; 



(* Expand necessary for later simplifications *) 



efil 
efi2 
efi3 
efi4 
efi5 



Expand [fil] 
Expand [fi2] 
Expand [fi3] 
Expand [fi4] 
Expand [fi5] 



efel 
efe2 
ef e3 
efe4 
ef e5 



Expand [fel] 
Expand [fe2] 
Expand [fe3] 
Expand [fe4] 
Expand [fe5] 
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A. 2 Mathematica Code for Linearizing Equations 
of Motion 

Following is the Mathematica code for analytically linearizing the non-linear equations 
of motion of the previous section. 

(* 

file for linearizing man5dof _sym_antisym model for following case. 

here the arm movement is decomposed into symmetric and 
antisymmetric components, phil = x4 + x5 , phir = x4 - x5 
x4 is the antisymmetric component 
x5 is the symmetric component 

vector from left limb e.g. to shoulder is {0,0, bj3} in limb coord. 

vector from right limb e.g. to shoulder is {0,0, bj3} in limb coord, 
mass of left and right limb are equal, mr = ml = m 
inertia of left and right limb are equal, Il=Ir 

vector from body e.g. to left shoulder is {0 , ibj2 , ibj3}in body coord 

vector from body e.g. to rt . shoulder is {0 ,-ibj2 , ibj3}in body coord 

nominal body rotation rate is {0,W,0} 

left arm is nominally at phiO, right arm is nominally at -phiO, 
deviations in both arm movements are of equal sign and magnitude. 

After the 5 dof linearized model is complete, a constraint equation 

of the form x = B z can be applied to the system of equations. 
B represents the dependence of the original state x on another 
state vector z. z can be of lower dimension than x. Then the 
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transformed set of linear equations are produced as follows: 

B~T M B zddot + B~T C B zdot + B~T K B z 

This transformation can be applied in this file or in a later 
one that evaluates the numerical quantities. 

*) 

(* for yman model *) 

(* 

simplify = {bjll->0, bjl2->0, bjl3->bj3,\ 

bjrl->0, bjr2->0, bjr3->bj3,\ 

ibjll->0, ibjl2-> ibj2, ibjl3->ibj3,\ 

ibjrl->0, ibjr2-> -ibj2, ibjr3->ibj3,\ 

irll->illl, ir22->il22, ir33->il33, mr->m, ml->m,\ 

mt->mb+2m} ; 
*) 
(* for yman model with no shoulder width and pvt at e.g.*) 

simplify = {bjll->0, bjl2->0, bjl3->bj3,\ 
bjrl->0, bjr2->0, bjr3->bj3,\ 
ibjll->0, ibjl2-> 0, ibjl3->0,\ 
ibjrl->0, ibjr2-> 0, ibjr3->0,\ 

irll->illl, ir22->il22, ir33->il33, mr->m, ml->m,\ 
mt->mb+2m} ; 



sfil = efil/ . simplify; 
sfi2 = efi2/ . simplify; 
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sfi3 = efi3/ . simplify; 
sfi4 = efi4/ . simplify; 
sfi5 = efi5/ . simplify; 

(* Replace angular velocities and accelerations with corresponding 

euler angle expressions. 
*) 

S213 = { 
{0 , Cos [thl [t] ] , Sin [thl [t] ] Cos [th2 [t] ] } , 
{0 , -Sin [thl [t] ] , Cos [thl [t] ] Cos [th2 [t] ] } , 
{1, 0, -Sin[th2[t]]} 

>; 

C213 = { 
{ Cos[thl[t]]Cos[th3[t]] + Sin[thl[t]]Sin[th2[t]]Sin[th3[t]] , 

Sin [thl [t] ] Cos [th2 [t] ] , 

-Cos[thl[t]]Sin[th3[t]] + Sin [thl [t] ] Sin[th2 [t] ] Cos [th3 [t] ] }, 
{ -Sin[thl[t]]Cos[th3[t]] + Cos [thl [t] ] Sin[th2 [t] ] Sin[th3 [t] ] , 

Cos[thl[t]]Cos[th2[t]] , 

Sin[thl[t]]Sin[th3[t]] + Cos [thl [t]] Sin[th2 [t]] Cos [th3[t]] } , 
{ Cos [th2 [t] ] Sin [th3 [t] ] , 

-Sin[th2[t]] , 

Cos [th2 [t] ] Cos [th3 [t] ] } 

>; 

thdot = {thl'Dt], th2'[t], th3'[t]}; 

wba = {0,W,0}; 

wca = S213. thdot + C213.wba 
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wcadot = D[wca,t]; 

(* See Hughs pp 116 and my notebood 6/14/93 pp 33 to see linearization 
of a steadily rotating system. This system steadily rotates 
about the 2 axis with speed W. 

the linearized states are not consistent with Hughs notation so 
that the mass matrix retains its nominal form. 

xl tilt angle (th2) 

x2 dev. from nominal sault angle (th3) 
x3 twist angle (thl) 

x4 dev of both arms in anit-symmetric mode form PhiO 
x5 dev of both arms in symmetric mode form PhiO 
*) 

eulerreplace = {wbl [t] ->wca[ [1] ] , wb2 [t] ->wca[ [2] ] , wb3 [t] ->wca[ [3] ] , 

wbl ' [t] ->wcadot [ [1] ] , wb2 ' [t] ->wcadot [ [2] ] , wb3 ' [t] ->wcadot [ [3] ] } ; 

linname = {phia[t] -> x4, phis[t] -> x5, \ 

thl[t] -> x3, th2[t] -> xl, th3[t] -> x2, \ 

thl'[t] -> xdot3 W, th2'[t] -> xdotl W, th3'[t] -> xdot2 W, \ 

phia' [t] -> xdot4 W, phis' [t] -> xdot5 W, \ 

thl"Et] -> xddot3 (W W), th2"Et] -> xddotl (W W) , th3"Et] -> xddot2 (WW), \ 

phia"Et] -> xddot4 (WW), \ 

phis"Et] -> xddot5 (W W)}; 

lincond = {xdotl->0, xdot2->0, xdot3->0, xdot4->0, xdot5->0,\ 
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xl->0, x2->0, x3->0, x4->0, x5->phi0}; 



eulfil = Expand [sf il/ . eulerreplace] 
eulfi2 = Expand [sfi2/ . eulerreplace] 
eulfi3 = Expand [sfi3/ . eulerreplace] 
eulfi4 = Expand [sfi4/ . eulerreplace] 
eulfi5 = Expand [sfi5/ . eulerreplace] 



renfil = Expand [eulf il/ .linname] 
renfi2 = Expand [eulf i2/ .linname] 
renfi3 = Expand [eulf i3/ .linname] 
renfi4 = Expand [eulf i4/ .linname] 
renfi5 = Expand [eulf i5/ .linname] 



renfel = Expand [ef el/ . linname] 
renfe2 = Expand [efe2/ . linname] 
renfe3 = Expand [efe3/ . linname] 
renfe4 = Expand [efe4/ . linname] 
renfe5 = Expand [efe5/ . linname] 



mmll = Coefficient [renf il ,xddotl] /. lincond; 

mml2 = Coefficient [renf il ,xddot2] /. lincond; 

mml3 = Coefficient [renf il ,xddot3] /. lincond; 

mml4 = Coefficient [renfil ,xddot4] /. lincond; 

mml5 = Coefficient [renfil ,xddot5] /. lincond; 



mm21 = Coefficient [renf i2,xddotl] /. lincond; 
mm22 = Coefficient [renf i2,xddot2] /. lincond; 
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mm23 = Coefficient [renf i2,xddot3] /. lincond; 
mm24 = Coefficient [renf i2,xddot4] /. lincond; 
mm25 = Coefficient [renf i2,xddot5] /. lincond; 

mm31 = Coefficient [renf i3,xddotl] /. lincond; 

mm32 = Coefficient [renf i3,xddot2] /. lincond; 

mm33 = Coefficient [renf i3,xddot3] /. lincond; 

mm34 = Coefficient [renf i3,xddot4] /. lincond; 

mm35 = Coefficient [renf i3,xddot5] /. lincond; 

mm41 = Coefficient [renf i4,xddotl] /. lincond; 

mm42 = Coefficient [renf i4,xddot2] /. lincond; 

mm43 = Coefficient [renf i4,xddot3] /. lincond; 

mm44 = Coefficient [renf i4,xddot4] /. lincond; 

mm45 = Coefficient [renf i4,xddot5] /. lincond; 

mm51 = Coefficient [renf i5,xddotl] /. lincond; 

mm52 = Coefficient [renf i5,xddot2] /. lincond; 

mm53 = Coefficient [renf i5,xddot3] /. lincond; 

mm54 = Coefficient [renf i5,xddot4] /. lincond; 

mm55 = Coefficient [renf i5,xddot5] /. lincond; 

(* this gets the d.c. forcing part *) 

del = renf il/ . lincond; 
dc2 = renf i2/ . lincond; 
dc3 = renf i3/ . lincond; 
dc4 = renf i4/ . lincond; 
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dc5 = renf i5/ . lincond; 



cl = renf il/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} 
c2 = renf ±2/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} 
c3 = renf ±3/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} 
c4 = renf i4/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} 
c5 = renf i5/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} 

dcdxll = D[cl ,xl] /.lincond; 

dcdxl2 = D[cl,x2]/. lincond; 

dcdxl3 = D[cl,x3]/. lincond; 

dcdxl4 = D[cl,x4]/. lincond; 

dcdxl5 = D[cl ,x5] /.lincond; 

dcdx21 = D[c2,xl]/. lincond; 
dcdx22 = D[c2,x2]/. lincond; 
dcdx23 = D [c2 ,x3] /. lincond; 
dcdx24 = D[c2,x4]/. lincond; 
dcdx25 = D [c2 ,x5] /. lincond; 

dcdx31 = D[c3,xl]/. lincond; 
dcdx32 = D [c3,x2] /. lincond; 
dcdx33 = D [c3,x3] /. lincond; 
dcdx34 = D[c3,x4]/. lincond; 
dcdx35 = D [c3,x5] /. lincond; 

dcdx41 = D[c4,xl]/. lincond; 
dcdx42 = D[c4,x2]/. lincond; 
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dcdx43 = D[c4,x3]/.lincond; 

dcdx44 = D[c4,x4]/. lincond; 

dcdx45 = D[c4,x5]/.lincond; 

dcdx51 = D[c5,xl] /.lincond; 

dcdx52 = D [c5 ,x2] /. lincond; 

dcdx53 = D [c5 ,x3] /. lincond; 

dcdx54 = D[c5,x4]/. lincond; 

dcdx55 = D [c5 ,x5] /. lincond; 



dcdxdotll = D[cl ,xdotl] /.lincond; 
dcdxdotl2 = D[cl ,xdot2] /.lincond; 
dcdxdotl3 = D [cl ,xdot3] /. lincond; 
dcdxdotl4 = D[cl,xdot4]/. lincond; 
dcdxdotl5 = D [cl ,xdot5] /. lincond; 

dcdxdot21 = D[c2,xdotl] /.lincond; 
dcdxdot22 = D [c2 ,xdot2] /. lincond; 
dcdxdot23 = D [c2 ,xdot3] /. lincond; 
dcdxdot24 = D [c2 ,xdot4] /. lincond; 
dcdxdot25 = D [c2 ,xdot5] /. lincond; 

dcdxdot31 = D [c3,xdotl] /. lincond; 
dcdxdot32 = D [c3,xdot2] /. lincond; 
dcdxdot33 = D [c3,xdot3] /. lincond; 
dcdxdot34 = D [c3,xdot4] /. lincond; 
dcdxdot35 = D [c3,xdot5] /. lincond; 
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dcdxdot41 = D[c4,xdotl]/.lincond; 

dcdxdot42 = D [c4,xdot2] /. lincond; 

dcdxdot43 = D [c4,xdot3] /. lincond; 

dcdxdot44 = D [c4,xdot4] /. lincond; 

dcdxdot45 = D [c4,xdot5] /. lincond; 

dcdxdot51 = D [c5 ,xdotl] /. lincond; 

dcdxdot52 = D [c5 ,xdot2] /. lincond; 

dcdxdot53 = D [c5 ,xdot3] /. lincond; 

dcdxdot54 = D [c5 ,xdot4] /. lincond; 

dcdxdot55 = D [c5 ,xdot5] /. lincond; 



(* this is the linearization of the external forces *) 

dtaudxdotll = D [renf el ,xdotl] / .lincond; 
dtaudxdotl2 = D [renf el ,xdot2] / .lincond; 
dtaudxdotl3 = D [renf el ,xdot3] / .lincond; 
dtaudxdotl4 = D [renf el ,xdot4] / .lincond; 
dtaudxdotl5 = D [renf el ,xdot5] / .lincond; 

dtaudxdot21 = D [renf e2 ,xdotl] / .lincond; 
dtaudxdot22 = D [renf e2 ,xdot2] / .lincond; 
dtaudxdot23 = D [renf e2 ,xdot3] / .lincond; 
dtaudxdot24 = D [renf e2 ,xdot4] / .lincond; 
dtaudxdot25 = D [renf e2 ,xdot5] / .lincond; 
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dtaudxdot31 = D [renf e3,xdotl] / .lincond; 
dtaudxdot32 = D [renf e3,xdot2] / .lincond; 
dtaudxdot33 = D [renf e3,xdot3] / .lincond; 
dtaudxdot34 = D [renf e3,xdot4] / .lincond; 
dtaudxdot35 = D [renf e3,xdot5] / .lincond; 

dtaudxdot41 = D [renf e4,xdotl] / .lincond; 
dtaudxdot42 = D [renf e4,xdot2] / .lincond; 
dtaudxdot43 = D [renf e4,xdot3] / .lincond; 
dtaudxdot44 = D [renf e4,xdot4] / .lincond; 
dtaudxdot45 = D [renf e4,xdot5] / .lincond; 

dtaudxdot51 = D [renf e5 ,xdotl] / .lincond; 
dtaudxdot52 = D [renf e5 ,xdot2] / .lincond; 
dtaudxdot53 = D [renf e5 ,xdot3] / .lincond; 
dtaudxdot54 = D [renf e5 ,xdot4] / .lincond; 
dtaudxdot55 = D [renf e5 ,xdot5] / .lincond; 



dtaudxll = D[renf el ,xl] /.lincond; 
dtaudxl2 = D [renf el ,x2] /. lincond; 
dtaudxl3 = D [renf el ,x3] /. lincond; 
dtaudxl4 = D [renf el ,x4] /. lincond; 
dtaudxl5 = D [renf el ,x5] /. lincond; 

dtaudx21 = D [renf e2 ,xl] /. lincond; 
dtaudx22 = D [renf e2 ,x2] /. lincond; 
dtaudx23 = D [renf e2 ,x3] /. lincond; 
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dtaudx24 = D [renf e2 ,x4] / . lincond; 
dtaudx25 = D [renf e2 ,x5] /. lincond; 

dtaudx31 = D [renf e3,xl] /. lincond; 
dtaudx32 = D [renf e3,x2] /. lincond; 
dtaudx33 = D [renf e3,x3] /. lincond; 
dtaudx34 = D [renf e3,x4] /. lincond; 
dtaudx35 = D [renf e3,x5] /. lincond; 

dtaudx41 = D [renf e4,xl] /. lincond; 
dtaudx42 = D [renf e4,x2] /. lincond; 
dtaudx43 = D [renf e4,x3] /. lincond; 
dtaudx44 = D [renf e4,x4] /. lincond; 
dtaudx45 = D [renf e4,x5] /. lincond; 

dtaudx51 = D [renf e5 ,xl] /. lincond; 
dtaudx52 = D [renf e5 ,x2] /. lincond; 
dtaudx53 = D [renf e5 ,x3] /. lincond; 
dtaudx54 = D [renf e5 ,x4] /. lincond; 
dtaudx55 = D [renf e5 ,x5] /. lincond; 
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A. 3 Parameters of the Five D.O.F. Linear Model 



In this section I present the algebraic expressions lor the elements ol the linear dy- 
namic model used to study passive layout stability. These elements were computed 
from the symbolic linearization of the full non-linear dynamic equations of the five 
d.o.f. system. The equations of motion and the associated matrices are repeated from 
Chapter 4 Equations 4.1-4.6. 

Mx + (G + D) x + ft 2 (K + K') x = 



where 



M 



G 



D 



mn 





"H3 











m 2 2 











™13 





?^33 

















m 44 

















m 55 





012 











-012 





023 











-023 























045 











-045 










2b' sh /tt 



2b' sh /tt 
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K 



&11 k 13 

k 22 

fei k 33 



& 55 



K' 







2k sh /n 2 





oo o o 2k sh /n 2 

The state vector, x, is comprised of (in order): f3 - the tilt angle, 7 - the twist angle, 
<f) a - the asymmetric deviation of the two arms from the nominal arm angle, a - the 
deviation of the somersault angle from the frame steadily rotating at rate 0, and <f> s 
- the symmetric deviation of the two arms from the nominal arm angle. 
The individual elements of the above matrices are given below. 



mn 



(Ib t + 2 Il t + 2 rb 2 2 mi + 2 rb 3 2 mi + 
8rb 3 2 mf 



8 rb 3 mf 4 rb 3 mf raj 



2 m; + raj 



4 r/ 5 rb 3 mi cos(<P 



(2m ; + raj) (2m ; + raj) 
16 rl 3 rb 3 mf cos(<P ) 
(2 m/ + raj) 



8 rl 3 rb 3 mf raj cos(<P ) 16 rl 3 rb 3 mf cos(<P ) 



(2 m/ + raj' 



+ 



2 m/ + raj 



+ 2r/ 5 2 m;cos(<2> ) 2 + 



8 rl 3 mf cos((Po) 4 r/ 5 m 2 raj cos(<P )^ 

\2 ' 



(2 m/ + m b) 



(2 mi + raj, 



8 rl 3 mf cos(([>o) . , = l2 • / .* \2^ 

+ 4 r/ 5 rb% mi sm( <Po ) + 2 r/ 5 m/ sm( <Po ) , 



2 m/ + raj 



mi3 = (2 77/ — 2 r/ 5 r& 5 m/ cos(<P ) 



8 r/ 5 r& 5 771}* cos(<P ) 



(2 m/ + raj) 
4 rl 3 rb 3 mf raj cos(<P ) 8 rl 3 rb 3 mf cos(<P ) 



(2 m/ + raj' 



+ 



2 mi + raj 



+ 2 r/ 5 2 m; cos(<2> ) 2 + 
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8rl 3 2 mf cos($ ) 2 4 rl 2 mf m b cos($ ) 2 8 rl 3 2 mf cos($ ) 2 

(2 m; + m t ) 2 (2 m; + raj) 2 2 m ; + raj 

2 r/ 5 r6 2 m; sin( <£>o ) + 2 r/ 5 m; sin( ^o) ) 

m 22 = (Ib 3 +2rb 2 2 mi + 2Il 3 cos(<Po) + 4 r/ 5 r& 2 ra; sin(# ) + 
2// 2 sin(^) 2 + 2r/ 5 2 m /S in(^) 2 ) 



m 33 = (2//, +2r/ 5 2 micos(<2> ) 2 + — ^ * L^- + 

(2 m/ + raj) 

Arl 3 2 mf m b cos(<P ) 2 8 rl 3 2 mf cos (<P ) 2 2 . 2 , 

;- ~2 7\ ; 1" 2 rl 3 m i sm ( $o ) , 

(2mi + m h y 2mi + m b 

o 8 r& 5 mf 4 r& 5 m? mj 8 rb 3 mf 

m 44 = (Ib 2 + 2 r6 5 m ; + -^ + 



(2 mi + m b ) (2 m; + mj) 2 m ; + mj 

16 r/ 5 r& 5 mf cos (# ) 8 r/ 5 rb 3 mf m h cos($ ) 

4 r/ 5 r& 5 m, cos(<2> ) ■ - 2 ■ - 2 + 

(2mi + m b ) (2mi + m b ) 

16 rU rbs mf cos((Pn) „ . , , 9 , 9 . , , 9 

' ^^ + 2 // 2 cos ( <2> ) 2 + 2 r/ 5 2 mi cos (<2> ) 2 + 

8 rl 3 2 mf cos($ ) 2 i rl 3 2 mf m b cos($ ) 2 8 rl 2 mf cos($ ) 2 

(2 m; + m t ) 2 (2 mi + m b ) 2 2 m ; + raj 

2//,sin(<2>„) 2 ) 



,9 , , 9 , 9 n9 8 rh mfsm(<Pn) 

m 55 = (2//, +2rl 3 2 m l co S (<P ) 2 + 2rl 3 2 m l sm(<P ) 2 + —^ '■ y —£L- + 

(2 m/ + mi) 

2_2 _ • /,fi \2 c , 2 9 



4 r/ 5 mf m b sin((Po) 8rl 3 mfsin{$oY 



(2mi + m b ) 2mi + m b 

gi2 = (Ibi - lb 2 + Ib 3 + 2 Il t +irb 2 2 m l -2 Il 2 cos( <P ) 2 + 

2 Il 3 cos($ ) 2 + 8 r/ 5 r& 2 m; sin( $ ) + 2 iZ a sin( $ f ~ 2 i? 5 sin( $ ) 2 + 
4 rl 3 mi sm(<P ) ) 
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5-23 = (-211! +2 Il 2 cos( <P ) 2 -2Il 3 cos( <P ) 2 -4 rl 3 rb 2 mi sm(<P ' 

2 , n r; • /a \2 , ; 2 • /* \2\ 



2// 2 sin($ ) +2// 5 sin($ ) - 4 r/ 5 2 m;sin($ 



. 16 r/ 5 rb s m^ sin (# ) 8 rl s rb s m 2 raj sin($ 
545 = (4 r/ 5 r& 5 mi sm( <2> ) + - -^ '- + : -^ ■ 



^n 



<43 



(2 m; + raj) (2 m; + raj 

16 rU rbs rn 2 sin(<P ) , , , , 

—^^ L^Z_ 4 // 2 cos ^ sin <P +4iZ 5 cos £ sin #o 

2 m; + raj 



2 . 16 r/ 5 2 mf cos($ ) sin($ ) 

4 rl s 2 mi cos(<2> ) sin(<Z> ) -i ^ ^— ^ '-- 

(2 mi + raj) 

8rl 3 2 mfm b cos($ ) sin(# ) 16 r/ 5 2 m 2 cos($ ) sin(# ) 

(2 m; + raj) 2 2 m; + raj 

,9 ,9 8 r& ? m? 4 rbs mf mi, 

(Ib 2 -Ib 3 -2 rb 2 2 mi + 2rb 3 2 mi + s - <-j + — - -i L* 

(2 m; + raj) (2 m; + raj) 

8rb 3 2 mf 16 r/ 5 rb s m^ cos($ ) 
- 4r/ 5 r& 5 m; cos($ ) — - j - 



2 m; + raj (2 m; + raj^ 

8 rl s rb s m 2 raj cos($ ) 16 r/ 5 r& 5 m 2 cos($ ) 2 

H ~ : r ^ «2 cos( <Po ) 



2 m; + raj) 2 m; + raj 



n9 , 9 , n9 8rf» m?cos(<P ) 

2 Il 3 cos(<2> ) 2 + 2 r/ 5 2 m;cos(<2> ) 2 + — - '- K -f- + 

(2 mi + raj) 

4r/ 5 2 m 2 rajcos($ ) 8 r/ 5 2 m 2 cos ( # ) . 

— ^ -^ — ^^ - 4 r/ 5 rb 2 mi sin (P 8 , 

(2 m; + raj) 2m; + raj 

2// 2 sin(<2> ) 2 + 2// 5 sin(<2> ) 2 - 2 r/ 5 2 m;sin(<2> ) 2 ) 



8rl 3 rbsm^ cos(<P ) 

(-2 r/ 5 r& 5 m, cos(<2> ) '■ -^ '- - 

(2 m; + raj) 

4rl 3 rb 3 m\ raj cos($ ) 8 r/ 5 r& 5 m 2 cos($ ) 2 

9 r r ^ "2 COS( Wg j — 

(2 m; + raj) 2m;+raj 

n9 , 9 , n9 8rf» mf cos(<Pn) 

2 // 5 cos ( <2> ) 2 + 2 r/ 5 2 m; cos (<2> ) 2 + — - * ^- + 

(2 m; + raj) 

4r/ 5 2 m 2 rajcos($ ) 8 r/ 5 2 m 2 cos ( # ) . 

— ^ -^ — '—^ - 2 rl s rb 2 mi sin ^ 

(2 m; + raj) 2m; + raj 

2// 2 sin(<2> ) 2 + 2// 5 sin(<2> ) 2 - 2 rl s 2 m;sin(<2> ) 2 ) 
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k 22 = {-Ibt + lb 2 -2 Il t - 2 rb 2 2 mi + 2 Ik cos( $ f - 4 rl s rb 2 m t sin( $ ) + 
2// 5 sin(<2> ) 2 -2rl 3 2 mism($ ) 2 ) 



8 r/ 5 r& 5 mf cos (# ) i rl s rb s mf m b cos($ ) 

k 33 = (-2 rl 3 rb 3 m, cos(<2> ) ■ -^ ■ -^ + 

(^mj + rnj) (2mi + m b ) 

\ 3 ' K —^- + 2 Il 2 cos {$ f -2 Il 3 cos($ ) + 2 rl 3 2 m,cos($ ) + 

2 mi + raj 

8 rl 3 2 mf cos($ ) 2 i rl 3 2 mf m b cos($ ) 2 8 rl 3 2 mf cos($ ) 2 



"^55 



+ 



(2 m; + raj) 2 (2 m ; + raj) 2 2 m ; + raj 

2iZ a sin(<2> ) 2 + 2// 5 sin(<2> ) 2 - 2 rl 2 m x sin(<2> ) 2 ) 



8rl 3 rb 3 mf cos((I> ) i rl s rb s mf m b cos($ ) 

(-2 r/ 5 r& 5 m, cos(<2> ) ■ - 2 ■ - 2 + 

(2mi + raj) (2mi + raj) 

8 rk rbs mf cos((Pn) , , ,9 , , ,9 , 9 , ,9 

* * ' K —^- + 2 Ik cos #o -2 Il 3 cos <2>o + 2 r/ 5 2 m,cos($ ) + 

2 mi + raj 

8 rl 3 2 mf cos($ ) 2 i rl 3 2 mf m b cos($ ) 2 8 rl 2 mf cos($ ) 2 



+ 



(2 m; + raj) 2 (2 m ; + raj) 2 2 m ; + raj 

2// 2 sm(^) 2 + 2// 3 sin(^) 2 -2r/ 3 2 m f sin(^) 2 - 8r / /32mfsm( ^ )2 

(2 m/ + raj) 

4 r/ 5 mf m b sm(<P ) 8 rl 3 mf sm(<P ) 



2 mi + raj) 2mi+ mb 
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A. 4 Analytic Solution 

Solving lor the torque-lree motion ol a rigid body tumbling in space is one ol the 
classic problems ol dynamics. Given the limited value ol numerical integration as 
a means ol studying dynamics problems before the advent ol the digital computer, 
much effort was applied to rinding an analytic solution to the governing equations ol 
motion. The analytic solution provides the angular velocity ol the body in a body 
hxed principle axes system and the attitude ol the principle axis system with respect 
to an inertial coordinate system [Hughes 86]. The contribution that these analytic 
solutions offer today is a concise description ol the minimal number ol states and 
non-dimensional parameters that govern rigid body rotation. The analytic solutions 
also give rise to elegant geometric interpretations ol the motion that helps provide 
some intuition ol this otherwise complex movement. 

I present the analytic solution lor a tri-inertial body. A tri-inertial body is one 
in which each ol the principal inertias is distinct. Without loss ol generality lor 
the rest ol this section I will insist that, Ii > I 2 > I3. The analytic solution for 
rotation ol a tri-inertial body uses the Jacobian elliptic lunctions. Elliptic lunctions 
are a generalization ol the circular lunctions, sine and cosine, that originated Irom 
an effort to compute the circumference ol an ellipse [Bowman 61]. The Euler angles 
describing the body attitude with respect to inertial space can also be written using 
elliptic lunctions. 

Two integrals ol motion are essential in the analytic solution for torque tree rigid 
body motion. These represent conservation ol angular momentum h and conservation 
ol kinetic energy T. 

h* = Ilul + P 2 ul + Itul (A.l) 

2T = I x u? x + W 2 + /3CU3 2 (A.2) 

For a hxed angular momentum, the kinetic energy ol rotation is bounded above and 
below by the kinetic energy ol two pure spin solutions. The highest possible energy 
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state corresponds to rotation about the minor principal axis and the lowest energy 
state corresonds to rotation about the major principal axis. A continuum of solutions 
of rotational motion of a rigid body exist, each corresponding to a different energy 
level, between these two bounds, 



h 2 h 2 

— < T < — 
2/i " " 2/3 



The qualitative nature of the solution is divided into two regions between these bounds 
and separated by the solution of pure spin about the middle principal axis. The 
following parameter is useful in distinquishing between the two qualitatively different 
solutions. 

I 



2T 



I has units of inertia. The bounds on I corresponding to the bounds on energy are 



h>I>h 



I = /1, I = 7 2 , and I = I 3 correspond to the pure spin solutions. 

The maximum amplitude of the body fixed angular velocities can be obtained 
directly from the two integrals of motion. 



LO\ r 



^Sr, 



<^2r, 



^2r, 



<^2r, 
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IIl(h 


-Is) 


h- 


-I 


Ih(h 


-Is) 


Il- 


-I 



Ih(Il - h 



h 
h 



I -h 

lUh - Is) 



1/2 



1/2 



1/2 



1/2 





(A.4) 


(h>I> h) 


(A.5) 


= h) 


(A.6) 


(h>I> h) 


(A.7) 
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Next, the analytic solution for the angular velocity is given for the case of h > I > I 2 



u>\ = s\U\ m dn{r\ k) 
u 2 = s 2 u 2m sn(T]k) 
lq 3 = s 3 io 3m cn(T; k) 



(A.8) 

(A.9) 

(A.10) 



and similarly for I 2 > I > I 3 



{h-h)(I-h 
Ihhh 

{h-h)(h-I) 
{h-h)(I-h)_ 



1/2 



(t ~ t ) 



1/2 



(A.11) 



(A.12) 



u>\ = Sicui m cn(r; k) 
u 2 = s 2 u 2m sn(T] k) 
u 3 = s 3 u 3m dn(T; k) 



(A.13) 
(A.14) 
(A.15) 



k 



[h ~ h)(h - I) 



IhI 2 I 3 



(h-h)(I-h) 

(h-h)(h-I) 



1/2 



(t ~ t ) 



1/2 



(A.16) 



(A.17) 



The coefficients Si,s 2 and s 3 are all ±1 such that Si-s 2 s 3 = — 1, allowing four com- 
binations of signs. The Jacobian elliptic functions, (sn } cn and dn), depend on the 
scaled time, r and the modulus k. When k— > sn, cn } dn— > sin, cos, 1. 

To complete the solution of the tumbling motion we require the attitude of the 
body as a function of time. Since angular momentum is conserved, the inertially 
fixed angular momentum vector provides the reference with which to measure body 
attitude. The components of angular momentum as measured in the body fixed axis 
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system is given by 



h! 




IlUl 


h 2 


= 


h^2 


h 3 




I3U3 



(A.18) 



Another expression for this vector can be obtained by projecting the inertially fixed 
angular momentum vector onto the body fixed axes using the matrix Cu (2.2). Since 
the orientation of the inertial coordinate system is arbitrary, we choose it to oriented 
such that the angular momentum vector is oriented intirely along the inertial y axis. 
With this assumption in mind, equating the components of these two expressions for 
angular momentum gives 



1 2^2 

h<^3 



-S@h 



Solutions for two of the Euler angles are available immediately as 



(A.19) 







sin 



tan" 



1 
-1 


h \ 

~IlUl~ 

. h^2 - 



(A.20) 
(A.21) 



The solultion for $ requires significantly more work. First we use equation (2.4) to 
substitute for Ui in (A.19). 

h(SyC & $ + CyQ) = S^Ceh 
I 2 (CyC&$ - S^Q) = C^Ceh 
h(S e $ + V) = -S e h 

Solving for $ from the first two of these equations yields 



$ 



Hhsl + hCl 



(A.22) 
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To proceed we need an expression for C\ and S\. Using equation A. 21, we obtain 



tan* 



Using trigonometric identities we obtain 



2 (/l^ ) 2 



(h^y 



sin 



2 $ 



[hui] 



cos 2 * 



{h^y + (i 2 u 2 y 
(h^y + (i 2 u 2 y 



(A.23) 
(A. 24) 



Using these in equation A. 22 we obtain 



$ = h 



hu>l + I 2 uj\ 

(h^y + (i 2 u 2 y 



Using equation 2.15 for the denominator and equation 2.16 for the numerator, $ can 
be written as 

$ = h 



2T - I 3 to 2 



h 2 -(I 3 u 3 ) 2 _ 

At this point we have to specialize the solution for one of the two qualitatively 
different cases. First we solve for the case I\ > I > I 2 . Substituting for u 3 from 
equation A. 10 

$ = h 



2T-I 3 slu>l m cn* 

h 2 - I 2 s 2 u 2 m cn 2 



When we use the identity sn 2 + en 2 = 1 we obtain, 



$ 



L\ + L 2 sn 2 T 
L 3 + L4 sn 2 T 



(A.25) 



where 



L 2 



h(2T - I 3 ul 



(A.26) 
(A.27) 
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L 3 = h 2 -I 2 u 2 3m (A.28) 

U = lM m (A.29) 



Integrating the equation above for <f> will provide the somersault angle as a function 
of time. This integral can be written in terms of the Legendre elliptic integrals. 

In a similar fashion we can solve for the somersault angle for the second case 
I 2 > I > I3. In this case substitute u 3 from equation A. 15 into A. 4 and use the 
following identity: dn = 1 — k 2 sn 2 . This results in the same form for $ with new 



coefficients. 



where 



L\ + L 2 sn 2 T 
L 3 + L4 sn 2 T 



U = h(2T-I 3 u 2 3m ) (A.30) 

L 2 = hI 3 to 2 3m k 2 (A.31) 

L 3 = h 2 -I 2 3 ul m (A.32) 

U = lW 3m k 2 (A.33) 



To integrate equation A. 25 we split the integrand into two parts to obtain 

L x f dt L 2 f^ nsn 2 dt 

$ = -r- / T~. ^ + 



L 3 J 1 + n sn 2 L 3 n Jto 1 + n sn 2 

where n = L\jL 3 . We would like to integrate this equation with respect to the scaled 
time r = C T t (A. 11 or A. 16) since sn = sn(r). Therefore, we compute the differential 
dt = tt- and substitute it into the integral 



L\ f dr L 2 /"*/ nsn 2 dr 

«= -r-^r / t-. ^ + 



L 3 C T J 1 + n sn 2 L 3 C T n Jto 1 + n sn 2 
To find the change in somersault angle in a set period of time we define the limits of 
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integration to be the initial time, t , and the final time, tf or equivalently r and tj. 
The change in somersault angle is then given by the following integral 



OLj — OLq 



L\ f T f dr L 2 f T f nsn 2 dT 



L 3 C T Jt 1 + n sn 2 L 3 C T n Jt 1 + n sn 2 



These integrals can be expressed as standard elliptic integrals 



dr 



and 



where sin < 



o 1 + n sn 2 (r) 



n sn 2 (r)dT 



n(&, n, <f>) 



t — n(&, n, 



/o 1 + n sti 2 {t) 
stit. Therefore, breaking up the limits of integration we obtain, 



OLf — CtQ 



L 3 C T 



dr 



o 1 + nsn 2 Jo 1 + nsn 2 



t <j dr 



u 



L 3 C T n 



f T f nsn 2 dr 



Jo 1 + nsn 2 



f T ° nsn 2 dr 



1 + ns 



:rr 



which is written compactly as, 



OLf-ao = -—— [Il(k, n, <f) f ) - Il(k, n, <^o)]+-j— ^ — [(t/ - H{k, n, <f> f )) - (r - Il(k, n, <f> ))] 



L 3 C T 



L 3 C T n 



To compute these integrals we must solve for the initial scaled time, r , for a given 
initial state. To this end, the angular velocity, cu 2 , is given by the elliptic function sn 

^2(7") = s 2 io 2m -sn(T; k) 



For a given value of cu 2 /(cu 2m s 2 ) we need to compute r. If we assume that r = 
when lo 2 is zero and increasing then by definition s 2 = 1. The inverse of the Jacobian 
elliptic function, sn, is the Jacobian elliptic integral of the first kind. 

>™ So UjJL 



^(1 - ;c 2 )(l - Pic 2 ) 
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Computation of this integral allows us to compute r for a given initial cu 2 . 
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